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Topological Electronic Liquids 

Electronic Physics of One Dimension Beyond the One Dimension 

Abstract 

There is a class of electronic liquids in dimensions greater than one, which show all essential 
properties of one dimensional electronic physics. These are topological liquids - correlated elec- 
tronic systems with a spectral flow. Compressible topological electronic liquids are superffuids. 

In this paper we present a study of a conventional model of a topological superfluid in two 
spatial dimensions. This model is thought to be relevant to a doped Mott insulator. We show 
how the spectral flow leads to the superfluid hydrodynamics and how the Orthogonality Catas- 
trophe affects off-diagonal matrix elements. We also compute the major electronic correlation 
functions. Among them arc the spectral function, the pair wave function and various tunneling 
amplitudes. To compute correlation functions we develop a method of current algebra - an 
extension of the bosonization technique of one spatial dimension. 

In order to emphasize a similarity between electronic liquids in one dimension and topological 
liquids in dimensions greater than one, we first review the Frohlich-Peierls mechanism of ideal 
conductivity in one dimension and then extend the physics and the methods into two spatial 
dimension. 

The paper is based on the lectures given at Troisieme cycle de la physique en Suisse romande 
EPFL, Lausanne June 1997. 
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1 Introduction 



Over last two decades the physics community has been and continues to be fascinated by the 
phenomena of one-dimensional electronic physics. Being so different from perturbative Fermi Liquid 
picture of metals, electronic physics in one dimension is robust and in many ways more consistent 
than the Fermi liquid itself. This tempts one to think that the phenomena we observe in one 
dimension are more universal than they appear now and must be found beyond the one dimension. 
This paper is intended to show that there is a class of electronic liquids in dimensions greater than 
one which exhibits all essential features of one dimensional physics. 

First, we must identify a core phenomenon responsible for the peculiar infrared properties of 
one dimensional electronic systems. We will argue that (i) this is the spectral flow (rather than a 
restrictive geometry) and moreover, (ii)that an electronic system in any dimension which exhibits 
the spectral flow, shows essentially the same property as a one dimensional system, and its low 
energy physics is essentially one dimensional. 

We refer to the electronic systems with a spectral flow as topological electronic liquids. Topo- 
logical electronic liquids are not exotic, in fact they are well known models of quantum field theory 
driven by an anomalous current algebra (examples are given below). In the last decade the evi- 
dence was mounting that topological liquids are relevant to strongly correlated electronic systems 
and may adequately describe the physics of a doped Mott insulator (another example of a physical 
system which may fall into a class of topological liquids that are quantum crystals around a melting 
transition). 

In this paper we will discuss compressible topological liquids. They possess an important prop- 
erty: their low energy sector is described by the hydrodynamics of an ideal quantum liquid. In 
dimensions greater than one, this alone means that compressible topological liquids are supercon- 
ductors. 

The spectral flow phenomenon leads to a number of anomalous features which make the physics 
of topological liquids different from the physics of the Fermi liquid and make the physics of the 
superconducting state different from BCS. The most dramatic is perhaps the Orthogonality Catas- 
trophe |l(|: an overlap between ground states differ by an odd number of particles tends to zero 
in a macroscopical system. This means that an electron is not a quasiparticle. While embedding 
an electron, the system rearranges its quantum state by creating soft modes with a singularly large 
density of states. As a consequence, most scattering amplitudes are singular at low momentum 
due to emission of a cloud of soft modes. This dramatically affects the physics of tunneling in the 
superconducting state. 

Here we discuss the the physics of the superconducting state of topological liquids and especially 
concentrating on various tunneling effects. 

On the technical side, the theoretical methods developed in one dimensional electronic physics 
may be extended to higher dimensional topological liquids. Among them is bosonization, a powerful 
tool used to compute correlation functions and matrix elements. We show how bosonization - a 
coherent state representation of a current algebra - works for higher dimensional topological fluids. 

Running ahead, let us give a simplified picture of the spectral flow phenomena fond found in 
compressible liquids. Let us consider gaped fermions in a static (vector or scalar) potential V. 
Assume that the chemical potential /i lies in a gap. When the potential V changes adiabatically, 
the fermionic energy levels shift. Generally the levels can not cross the energy E = [i. However, 
there are potentials whose adiabatic and arbitrary smooth variation may create some (always even) 
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number of unoccupied state below E = fi. or force some occupied levels to cross the level E = /i. 
This is the spectral flow or level crossing. To produce a level crossing, a variation of the potential 
( a soliton) must necessarily be topological. The index theorem then relates the topological charge 
of the potential and the number of levels crossed. If this simple phenomenon occurs, we conclude 
that a compressible liquid is a superconductor. Here is a short hierarchical list of familiar models 
that yield topological liquids: 

One Dimension. The Peierls model: 

H = i/)l(a x id x + /?7ri + ia x f5-K 2 )ipa (1-1) 

(7=1,2 

where the Dirac matrices a x ,f3 may be chosen as the Pauli matrices: a x = 03, /3 = o~\ and the 
modulus of the vector (7ri,7r2) is assumed to take a fixed value at infinity. 

The soliton here is a kink of the 7r(x)- field. 

Two Dimensions. Dirac fermions coupled to a vector field field 

H = ift(aiV + (3ttt + f3m)ip . (1.2) 
where r are Pauli matrices. Solitons in this model are hedgehogs (skyrmions) of the vector field ir. 

Three Dimensions. Dirac fermions coupled to vector chiral bosons: 

H = ^ (iaV + +i-f 5 (37TT + Pm)i; . (1.3) 

Solitons are skyrmions of the vector bose field. The last model is known as a linear a model of 
current algebra and describes nuclear forces. 

All these models in dimensions greater than one, being treated at fixed chemical potential, are 
superconductors! 

Below we will concentrate on another model of Dirac fermions interacting with a gauge field 
(not an electromagnetic field): 

H= ^(«HV + A) + /3m)^ (7 . (1.4) 

(7=1,2 

Although the physics of all the models is similar, this particular model is selected because of 
historical reasons and because one may trace it genesis from a model of doped Mott insulator. 
Solitons in this model are vortices of the gauge field. 

The text is organized as follows. We first review some basic elements of superconductivity with 
an emphasis on general properties which do not depend on the mechanism (Sec. 2). Then we review 
the Frohlich's ideal conductivity of the one dimensional model Ql.l| ) - a one dimensional topological 
liquid (Sec. 3). On this example we illustrate how the spectral flow leads to the hydro dynamics, and 
how the orthogonality catastrophe affects matrix elements and determines the physics at the Fermi 
surface. We also review methods of computing a mass shell asymptotic of correlation functions 
(bosonization) and develop a vertex operator technique to compute matrix elements. This, except 
the matrix elements, is a standard material. We adopted it for the purpose of generalization to 
higher spatial dimensions. 
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In the next section (Sec. 4), we extend the methods of Sec. 3 to the model ( |1.4[ ) of topological 
liquid of two spatial dimension. We try to follow the line of Sec. 3 as closely as possible, in order to 
separate the physics of one dimension which is caused by to a restrictive one dimensional geometry 
and the general physical aspects of spectral flow. 

Finally in Sec. 5 we discuss various tunneling mechanisms, where a signature of topologicall 
mechanism of superconductivity is the most transparent. 

A comment on literature is in order. We shall not give too many references for the one dimen- 
sional part of this text. The Frohlich's conductivity goes back to his original seminal paper of 1955 



[14], which impresses even todays readers. Excellent reviews on the further development of physics 
of charge density waves are fi~5j| . The part regarding the bosonization and correlation functions has 
been developed in late 70'th and redeveloped in early 90'th. Many classical papers may be found 
in an excellent preprint collection §. As the more recent paper we refer to ||. The vast body of 
Sec. 3 is in fact the PhD thesis of the author of 1980. The two-particle matrix element calculation 
of Sec. 3. 5 is believed to be new. Regarding the spectral flow, current algebra, anomalies and the 
geometrical phase there is also excellent preprint collections ||. Most of the material of the Sees. 4 
and 5 are based on recent papers of A. Abanov and the author ||, [|. We do not refer to recent 
attempts to extend the bosonization to higher dimension as an alternative description of the Fermi 
liquid, since they do not seem relevant to the subject of Sec.4. 

2 Criterion for Superconductivity. 

We start from a discussion of the definition, criteria, and implementation of superconductivity. 
This section is aimed to separate fundamental aspects of the phenomenon of superconductivity 
from features inherited through a particular(BCS) mechanism. 

2.1 Hydrodynamics of an Ideal Liquid Phenomenology of Superconductors 
2.1.1 Euler equations for an ideal liquid 

The phenomenology of superconductivity starts by identifying the superfluid with an ideal quantum 
liquid. The hydrodynamics of an ideal liquid are given by two equations M — continuity 

d tP + Vj = (2.1) 

and conservation of momentum (Euler equation) 

8 t ji + ^-U ik = 0, (2.2) 
dr k 

where 

= p5 ik + pviV k (2.3) 
is the momentum tensor, and p, v, j = pv,p are the density, velocity, current and pressure. 
All together they can be assembled in the Hamiltonian form 

B = (2.4) 

[p(r),v(r')] = -±V r <S(r-r') (2.5) 
m 
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where e{p} is the potential energy of the liquid. In the linear approximation p — p s <C p s (a 
super fluid density) and 



e M = |"P 2 , P = vl{p-ps), (2.6) 

where w is a compressibility and = dp/ dp = is the sound velocity, we obtain the Hamiltonian 
of linear hydrodynamics [] 

H = f K 2 j 2 + (p - Ps) 2 } + A cxt j + A^\p - Ps ) 

[p(r), j(r')] = -ip s V r S(r - r') (2.7) 

Here and in what follows we use k = wvq A = de/dp and also refer to it as a compressibility. We add 
an external electromagnetic potential in order to test a response of the system to an electromagnetic 
field. These give the basic equation of the superfluid — a quantum metastable state with a current: 



1 



Vxj = "4^ B ' (2 ' 8) 

L 

v ^ Vp + dt j = -L^B cxt . (2.9) 

The first equation (together with the Maxwell equation VB ext = 47rj) implies the Meissner 
effect, where 

A i = (4vreW- 2 )- 1/2 = (£,)V 2 (^) = cK 

is the London penetration depth (for London's superconductors) and u> p = (4ir p s e 2 / m) 2 is the 
plasma frequency. The second equation means that the liquid is compressible and has an ideal 
conductivity (a(uj) ~ iu;^ 1 ) 0. 

In the Lorentz gauge VA + VQ 2 dtAo = 0, these equations acquire familiar form 

' A cxt . (2.10) 



4vrAl 



The London equations can be also translated into the current-current correlators at small u 
and k. Varying these equations (2.£,2~!9l) over A ext and taking into account the continuity equation 

O we g et 

(p(q,u)p(- q -u;)) = -K- 1 fi q \ - (2.11) 



UJ 2 — VqO 2 



and 

-1„.2 



CLHZ, -w)j±(?, -w)) = k- l v$, Vj ± = , (2.12) 

where jj_ is the transverse component of the current. This is another way to say that the superfluid is 
a compressible ideal liquid. The density-density correlator ( 2.11Q shows a gapless mode (longitudinal 



sound) which means that the liquid is an ideal and compressible, i.e. may flow without dissipation. 

Parameters m and p a being treated separately are fictitious. Only their ratio, i.e. compressibility w = -J 5 — and 



the sound velocity Vq are physical and measurable. 

2 The Coulomb interaction between charge carriers makes the liquid incompressible. Below we neglect the Coulomb 
interaction and consider the electromagnetic field as external. 
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The second equation (which equivalently reflects the Meissner effect) implies that there is no gapless 
transverse mode, i.e. that there is no shear modulus. A zero rigidity to shear is the only true 
definition of a liquid. As a contrary transverse sound is a feature of a solid. For a solid the density- 
density correlator is the same is in eq. ( 2.1 1| ) , whereas for the transversal current correlation one 
would have 



2 



(j±(g,k>)j±(-g,-k>)) = -« ^0-5 — 2~2 ■ ( 2 - 13 ) 

In one dimension, there is no true distinction between solid and liquid — there is no room to apply 
a shear. 



2.1.2 Displacements and Potentials 

Density and current are not cannonical variables. In order to quantize the hydrodynamics canoni- 
cally, one need to solve continuity equation (|2.lD . The latter plays a role of a constraint. It can be 
done through displacement u: 

p-p s = -Vu, j = d t u, (2.14) 
which is a cannonical partner of the current 

fo(r), ufc(r')] = -A ik 5(r - r'). (2.15) 
m 

In the linear approximation the displacement u is a free field. In terms of displacements the linear 
hydrodynamics (2/7) becomes 

£ = ^(v 2 (d t uf-(Vuf). (2.16) 

The most familiar origin of superfluidity in quantum systems is a formation of a condensate. In 
this case the continuity equation can be also solved by setting 

j = 2^*(^V " 2eA ext )^ + h.c, 

[¥*(r),¥(0)]=<5(r). (2.17) 

where ^> = l^e 1 ^ is an amplitude of the condensate - a sort of order parameter. Then the hy- 
drodynamics can be written as Gross-Pitaevsky equations qI Substituting ( 2.1 7| ) into (2.6), we 
get 

id t m = —V 2 ^ + k{\^\ 2 - Ps )^ (2.18) 

2m* 

Parameter m* establishes the second scale for the superconductor - the correlation length £ ~ 
h/m^vo. 

If one neglects any inhomogeneity of the modulus \^/\ , the current becomes a gradient of the 
phase 

j = --^-V0. (2.19) 
2m* 



3 The Gross-Pitaevsky equation resembles, but has a little in common with the Ginzburg-Landau equation. The 
latter is valid near T c . The applicability of the former equation to superconductors, if any, is limited. 
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If for some reasons one may neglect vorticity (i.e. transversal part of the current), this representation 



(regardless the equation ( |2.18D ) gives the standard version form of the Josephson effect. It can be 
seen in the form of point current-current correlator 

(j(w)j(-u;)>~%). (2.20) 

This result does not follow from the hydrodynamics alone, but is rather based on a strong assump- 
tion that the phase <fi has no vorticity and depends on the mechanism. It different in topological 
superconductors. 

In fact vorticity is involved in a general solution of the continuity equation. Let us combine a 
density and a current to d+ 1 vector = (po,j). The continuity equation then requires the current 
to be a 1-form. In dimensions 1, 2, 3 it is respectively 

3(i = d(j,<i> + Z(xv\d v A x , 

3(i = d^4> + e^xpduFxp, (2.21) 

where F\ p , A\,<p are potentials. In one dimensions the potential tp is also a displacement. 

In contrast to the BCS, in topological liquids, potentials are more natural (less singular) fields 
than BCS phase eft. This obviously weakens the Josephson effect but does not eliminate it and the 
superconductivity - the current- current correlator ( 2.20| ) has a power law peak around co = with 



the width of the order of Fermi energy, rather than being a delta function. 
2.2 Landau Criterion 

Not every electronic liquid is ideal. A Fermi liquid, for example, are dissipative. Their current 
correlation functions show Landau damping: 

f k<k f m 

^m(j ± (q,uj)j ± (q,uj)) r / S(u + e(k) - e(k + q))dk ~ -u mkt- , (2.22) 

J\q.+k\>k f q 

where vq is a density of states. The origins of the dissipation are 

(i) gapless particle- hole excitations, and 

(ii) leaving the Fermi surface, an electron can easily change the direction of the momentum (the 



dissipation is due to the integration in ( 2.22j ) over the angle between momenta k and k + q). For 



obvious reasons there is no Landau dissipation in the one dimensional case. 
To avoid the dissipation a liquid either 

{\)must have no gapless density modes other than longitudinal sound. In particular the single 
particle spectrum of the homogeneous liquid must have a gap., and/or 

(ii) all scattering channels which change a direction of an electron must be effectively suppressed 



The first proposition is the subject of the Landau criterion Q. In its grotesque form it states 
that 

the spectrum must contain longitudinal sound and the single particle spectrum must have 
a gap. 
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In the mechanism we discuss below both these sufficient criteria hold true. 



All of this is true for a homogeneous liquid. The true check whether a liquid is superconductive 
may be made only if some non-zero concentration of impurities does not result in resistivity. There 
is a belief that if the Meissner effect holds in the absence of impurities, a weak disorder does not 
lead to resistivity, regardless of the mechanism of superconductivity or scattering. 

The Landau criterion seems to be sufficient in spatial dimensions greater than one. In one 
dimension it fails for a very simple reason, namely, in one dimension we cannot distinguish between 
a liquid and a solid since there is no shear and no Meissner effect. As a result a single impurity 
pins down the flow in the same way as a single impurity pins down a solid [JjJ. 

An important consequence of a gaped one particle spectrum and a soft density mode is that the 
gap opens always at the level of the chemical potential. This means that while adding any even 
number of particles into the system, the spectrum rearranges itself so that new particles appear 
below the gap. Then the two particle wave function, i.e., the matrix element 

A(n,r 2 ) = <AA + 2| C {( ri ) C {(r 2 )|A0 (2.23) 

between the ground states of the system with N and N + 2 particles is localized, i.e. decays with 
a distance between particles. This gives rise to the Josephson tunneling. 

The Josephson tunneling amplitude, i.e. the two particle matrix element must not be confused 
with the order parameter of See p. 1.2 . 



2.3 Implementation of superconductivity. 

Superconductivity — a metastable quantum state with a current in a macroscopical system — 
manifests itself as a particular set of correlations in the ground state: 

• (i) Meissner effect: 

no rigidity to shear (an ideal diamagnetism) 

(j±(k,uj)i ± (-k,-cj)) = (4 Js 2 , &tk,u^0. (2.24) 

• (ii) Gap in the electronic spectrum: a singularity uj = f2(p) closest to the origin of the one- 
particle Green's function in the u> plane 

G(w,p) = (4(cu,p) C(T (cu,p)) (2.25) 

• (iii) Josephson tunneling: the two particle wave function (the matrix element between the 
ground states of the system with N and N + 2 particles) is localized and normalizable 

A(n -r 2 ) = {N + 2\c\{v l )c\{v 2 )\N) (2.26) 

• (iv) Correlation length: 



(4vrA L ) 



;i - o(kzf)A 



cxt 



(2.27) 
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These correlations describe very different aspects of the phenomenon: (i, iv)the hydrodynamics 
of an ideal liquid (density and compressibility), (ii) one particle spectrum, and (hi) two particle 
matrix element, (iv) correction to the hydrodynamics. 

They also set up scales: a penetration depth, correlation length, a tunneling amplitude, a gap 
and a transition temperature T c . 

From a general point of view, there is almost no reasons to relate different quantities and 
different scales to each other. 

Nevertheless, due to the mean field character of the BCS theory many of them turn out to be 
essentially the same. For instance, the gap ( 2.25| ) and two particle matrix element ( p. 23 ) appear to 



be related and the two particle matrix element is identified with the order parameter ^. Moreover, 
the inverse correlation length, gap, tunneling amplitude and T c emerge as the same scale. 

This misleading "advantage" of the BCS theory often allows one to draw conclusions about the 
gap function by looking at the matrix element, about the tunneling based on the hydrodynamics, 
and about transition temperature based on the gap and vice versa. 

However, the gap, Josephson current, and the order parameter are essentially different: the first 
characterizes the spectrum, the second is a matrix element, determined also by the phase of the 
wave function, while the third measures correlation between pairs. For instance, in the model we 
consider below, contrary to the BCS, the pair correlation function 

lim <4(ra) C {(r 2 ) CT (r 3 )q(r 4 )> - (2.28) 

|(ri,r2)-(r3,r4)| >oo 

vanishes due to strong fluctuations of the phase of a pair et(ri)ct (r 2 ). This however does not mean 
that there is not a sort of a long range order or superfluid density ( 2.27| ) . 



These quantities start to diverge beyond the mean field level of the BCS theory, although the 
deviations are perturbative (non singular). 

In an electronic liquid where the interaction is strong, one also expects to see a difference between 
dissimilar implementations of superconductivity. This difference becomes dramatic in topological 
liquids, where the entire effect of superconductivity is due to peculiar quantum phases of wave 
functions of the ground state and low energy excitations. The deviations are also dramatic in the 
cuprate superconductors. 

There are several reasons why topological liquids are interesting. First of all they are new 
electronic liquids, known previously only in one dimension. In higher spatial dimensions they 
exhibit a fundamentally new mechanism of superconductivity and superfluidity. Secondly, the 
topological mechanism has been found in models of strongly correlated electronic systems, namely 
in the doped Mott insulator. 

Below we provide a detailed analysis of a model of topological superconductivity. 

We start by reviewing the one dimensional case were the topological character of the ground 
state is the most transparent. 

3 Prohlich Ideal Conductivity — One Dimension. 

In the early days of superconductivity and before BCS, Frohlich [Q noticed that in a one- 
dimensional metal an incommensurate charge density wave (CDW) slides through the lattice unat- 
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tenuated. Since it carries an electric charge and since a gap has been developed in the electronic 
spectrum, Frohlich concluded that the ground state of his system is superconducting. Although a 
sliding charge density wave indeed contributes to conductivity, the possibility of Frohlich supercon- 
ductivity in one dimension was considered nothing more than a theoretical curiosity, because of a 
variety of pinning mechanisms Q . 

The failure of the charge density wave mechanism in one dimension does not devaluate Frohlich's 
ideas, which as we shall see are valid in higher dimensions |D| |. 

Moreover, Peierls-Frohlich model reveals the physics of the bosonization - the method which 
work just as well in higher dimensional topological liquids. 



3.1 Peierls Instability 

Let us now review Frohlich's ideas (see e.g.,Ref. fl5|]). We start from the Frohlich model of an 
incommensurate electron-phonon system 

n ( V 2 \ 

H = \^ ( cjl( \i)c a + gu(x)c* V c a ) + an energy of phonons (3-1) 

, V 2m / 

<T=1 

where c a are electrons, u(x) is a displacement (phonons) and n is the degeneracy of the electronic 
state (spin, for example). 

In one-dimensional electron-phonon systems the Peierls instability causes a lattice displacement 

(u(x)) = — cos (2k f x + if) (3.2) 
9 

with an amplitude Ao ~ Et exp(— const /g 2 ) <C Ef and a period equal to the average distance 
between particles 2ir/2kf. In its turn the phonon modulation causes a modulation of the electronic 
density p — po ~ Ao cos(2/cpx + if). The phase <p determines the position of the charge density wave 
(CDW) relative to the lattice. 

Frohlich noticed that the periodic density fluctuations of electrons are fixed only relative to 
the lattice (clearly the energy of the incommensurate state does not depend on the constant part 
of (p) and it can easily travel with some velocity, such that p ~ Aqcos(2/cf(x — vt) + (p). This 
could be compensated by changing the phase according to <p = —2kpv. Therefore the current 
j = pv = v^nkF/n), is 



and from continuity 



P = A) + ^V • (3-4) 



These are Frohlich's equations. 

If the charge density wave is not pinned it propagates with the dispersion oj = v$k (where vq 
needs to be computed) i.e., 

£c = -^—((d t <p) 2 -v 2 (d x <p) 2 ), (3.5) 

where Vf = ^ / v$ is the Fermi velocity, and due to the Frohlich equations yields the hydrody- 
namics ( [2. 16 ). 
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Since there is a gap in the electronic spectrum, and the only gapless mode is the sliding CDW, 
Frohlich concluded that his system is superconductive. In fact the CDW is an ideal conductor 
rather than a superconductor due to pinning mechanisms. The Landau criterion does not work in 
one dimension - there is no space to flow around an obstacle. 

Let us now derive the Frohlich equations formally. The first step is to pick out the fast variables 
and keep only the slow variables. In this case the slow variables are associated with electrons in 
the vicinity of two Fermi points ±kp, 

c(x) ~ e ikFX ip L + e~ ikFX ij R , (3.6) 

and phonons with momentum close to 2k t 

A(x) 

u(x) ~ — — cos(2/c F x + (fix)) . (3.7) 
9 

where (p(x) and A(x) are fluctuating fields. In the continuum limit we then obtain the so-called 
linear cr-model 

• 2 2 

L = LL - Li- + ${ib - | A|e^)V > ( 3 - 8 ) 
g z oo z g z 

where D = %(i5 M — A^ xt ) and 70,71,75 are two-dimensional Dirac matrices, and uJ is a characteristic 
frequency of phonons. The modulus of the phonon field does not fluctuate much and is determined 
by its mean field value Ao- The effective model is 

L = ^ CT (iD-A e^)Vv , (3.9) 

This is the result of the Peierls instability - a gap Ao has opened at the Fermi level. What is 
the spectrum of this system? We shall see that the spectrum consists of 

• a gapped electronic mode 



e{k) =±y/A 2 + v 2 f (k±k f ) 2 ; 

• a gapless sonic mode - a phase ip of distortion. Below we shall see that it is a mode of 
modulation of density with dispersion E(k) = vok; 

• a gapped mode of spin density modulation. 

The vacuum of the model is infinitely degenerate: the states \<p) and \<p + const) have the same 
energy. Degeneracy leads to a soft sonic mode so that the system is an ideal conductor. Let us 
stress that position of the gap is always at the Fermi level, so that the spectrum strongly depends 
on the number of particles (filling factor). 

For comparison let us consider a commensurate Peierls model, where the number of electrons 
(with a given spin) is half the number of lattice sites. A canonical tight-binding model is 

H = J2 A n ,„+i(a+a n+ i + h.c) + H A (3.10) 

n 

where A n>ri+ i is a fluctuating hopping amplitude and is a phonon energy. In the continuum 
limit the half-filled case is described by the same a- model ( |3.8| ) but with a real A: 

: ,2 



L = ^(iD-i l5 A(x))^-^ 1 + ^ (3.11) 

2g z g ^ 
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At exact half filling, the CDW is two- fold commensurate and so the vacuum is two-fold degener- 
ate. In this case there is no soft translational mode since the CDW is commensurate and is pinned 
by the lattice. The excitations are gapped electrons and kinks of A (x) which connect two-fold 
degenerate mean field vacua: A — ► ±Ao when x — > ±00. The system is an insulator. It becomes 
an ideal conductor under a doping. 

■ 2 

Below we discard the kinetic energy of phonons (the term p=^- in eqs. (|3.§|,3.11)) by sending 
g — > 00. The energy of phonons then appears as radiative corrections due to exchange by electrons. 
The advantage of this limit is that the compressibility of the system (i.e. its central charge) and 
conformal dimensions of all operators will be universal and determined by the number of fermionic 
flavors n. Later we specify that n is equal 2. 



3.2 Hydrodynamics 

3.2.1 Spectral flow and Zero modes 

Let us now add particles to the system. In an incommensurate case the gap must follow the new 
chemical potential. In other words extra particles will not go into the upper band to occupy the 
lowest empty state. Instead a doping rearranges the period of the CDW, so as to create one more 
level in the lower band. This phenomenon is known as it level crossing or a spectral flow. How 
does it happen? 

In the presence of a kink in the phase, 92(00) — ip(-oo) = 2tt, the electronic spectrum changes in 
such a way that one extra unoccupied level appears at the top of the occupied band with an energy 
E = — Ao- When we add a particle to the system, it will therefore create a kink in the spatial 
configuration of p and an extra level, in order to be absorbed by the lower band. Moreover, the 
density of extra particles is locally and adiabatically bound to the kink as is suggested by eqs.( |3.3| , 

H). 

Tl 

jfi = -^e^duip . (3.12) 

We derive this equation on a more formal basis later on. 

The situation in the commensurate case ( |3.11|) is more subtle but is essentially the same. In 
the presence of the kink the electronic spectrum remains approximately unchanged, except for the 
appearance of the so-called zero mode, a state with a zero energy, located exactly in the middle of 
the gap. Indeed, the solution of the Schroedinger equation at a static kink A(x — > ±00) — > ±|Ao| 

id x ip L - iAi; R = E^ L (3.13) 

- id x ip R + iAtf; L = E^p R (3.14) 
with minimal energy is the zero mode 

ip L = ipR = e-Io Aix) ; E = 0. (3.15) 

The wave function of the zero mode is located in the core of the kink. 

While interpreting this result a subtle difference occurs, according to whether we keep the 
number of particles or the chemical potential fixed. If we fix the number of particles, the zero 
mode is unoccupied and ready to accept an extra particle. In the case of a fixed chemical potential 
(equal zero), the zero mode is occupied by 1/2. To understand this better let us notice that a single 
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kink (an odd number of kinks) is not compatible with fixed boundary conditions. If the boundary 
conditions are fixed, say periodic, then the minimum of two kinks (an even number of kinks) is 
required. However the two kinks may be well separated and can be treated almost independently. 
Since the topological charge is zero, there is no zero modes, but rather a symmetrical and an 
antisymmetrical combination of zero modes of independent kinks, split around zero energy. The 
antisymmetric state will appear below the chemical potential. If the chemical potential is kept at 
zero, it will be occupied by a particle. Thus we have two kinks per particle, i.e. each kink will be 
1/2 occupied. 

An adiabatic relation between the density of particles and soliton configuration (the axial current 
anomaly), similar to eq.(3.4), in this case tells us that the density of extra states is equal to the 



density of zeros of A(x) 

p(x) = ^(A&n))^ • (3.16) 

The factor 1/2 reflects the fractional occupation number, discussed above. 

Suppose we now dope the system by adding extra particles. The remarkable fact is that the 
system will lower its energy by creating the number of solitons (zeros in A(cc)) which is necessary 
to absorb all dopants. We refer this phenomenon as a topological instability. 

If a nonzero density of particles is added, i.e. a nonzero density of kinks is created, an interaction 
between kinks results in the formation of a narrow band in the middle of the wide gap (a midgap 
band). The width of this band is of the order of Ao exp(— const ^ ) where p — p s is the doping 
density and p s the density of the undoped system. This band absorbs all the dopants and is always 
completely filled. 

The passage to an incommensurate case is as follows. As in Prohlich's case, solitons have a 
translation mode due to their topological origin: a soliton lattice can slide along the atomic lattice 
without dissipation. Let X{ be the zeros of A(x). Then the density of extra particles (dopants) 
is Sp(x) = \ ^ZiS(x — Xi). Displacement of the positions of zeros around their mean field values, 
Xi = xi + ip(xi,t)/2irpQ, give rise to fluctuations of the density 5p(x) = p(x) — p$. According to 
( 3.16| ) they obey the same Frohlich equations ([T^, 3.4, 3.12| ). 



Each twist of <p adds one additional state in the middle of the gap. Therefore, adding n e extra 
particles gives rise to the topological charge of n e = Q where Q = f All of this is true in 
incommensurate cases when the system, after doping, has infinitely degenerate classical vacua. If 
the doping is a rational number, say, p/q, then the number of degenerate vacua is finite, namely q. 
The CDW is generally pinned by an exponentially small potential Aq(Ao/£/) 9-2 cos qcf). 

If the CDW is not pinned, an electric field can easily drag the solitons relative to the crystal. 
This leads to a current and is given by axial current anomaly equations, which determine the 
response to an external electromagnetic field. Setting m = Vq = 1, it reads 



e»M» = ^E cxt . (3.17) 



These equations are equivalent to the equations of linear hydrodynamics ( |2.8j , 2.9). As we 
discussed earlier, we cannot distinguish between a solid and a liquid in one dimension. Nevertheless, 
there is a global version of the Meissner effect (a spectral flow) in one-dimension- a diamagnetic 
current is generated by a magnetic flux $ ext being set inside a metallic ring: 



71 

jdx = $ ext . (3.18) 

7T 
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Combining ( |3.17| ) with the Frohlich equations ( 3.12| ) and using the relation j = —dL/dA cxt we 



obtain a bosonized version of the incommensurate CDW: 



L * = £{hM 2 -*0*<p) ■ ( 3 - 19 ) 



7T 



The Hamiltonian of the linear hydrodynamics (^?j) in one dimensional literature is known as 
Sugawara form. 

3.2.2 Bosonization 



For illustrative purposes let us rederive the results of the Sec, |3.2.1 by the bosonization procedure. 



For simplicity let us consider only two fermionic species (spin) a =|, j. In this approach, fermions 
ip a are treated as soliton of a boson field 4> a : 

^La ~ a- 1 ' 2 : e i4>L ° :, Vftx ~ a~ 1/2 : e - ^*" : (3.20) 

where : . . . : is the normal ordering and 4>l,r are holomorphic components of the cannonical Bose 
field 

[4>l{x), 4>L{y)] = - [4>r(x), <Pr(0)] = ivrsign(x - y), 

M%)MQ)- ■ M*)M0) ■ = In(-), (3.21) 

X 

where a is a lattice scale and L is the size of the system. The chiral components of currents in 
terms of chiral bosons are 

+ 1 +1 

]L<r = TpLa^La = ^d x (j) La , j R(T = ^ Ra ^Ra = -^dx4>Ra (3-22) 

whereas the Hamiltonian of free fermions is given by the Sugawara form: 

H = i^dxipLa - iiplfrdxipito = — ((d x <t> La ) 2 + (d^Rr) 2 ) . (3.23) 

It is instructive to rewrite these formulas in the Hamiltonian formalism, in terms of displacement 
and its cannonical momentum 

(pa = ^{4>L<7 + <t>Ra), R(T = \%{4>ha ~ <pRa), 

[0 CT (x),n CT (o)] =iS(x), 

p a {x) = d x cj) a {x), j a {x) = in,(x). (3.24) 
In these terms the electronic operator ( |3.6[) is 



C(T (x) ~ e »Arg(kf) e ik/a: . g i n CT (a/)dz'+i7rsign(k/)<£ CT (a;) . (3.25) 
k f 

where kj = ±kf and the amplitude of backward scattering is 

A# 75¥ > ~ A : cos(2tt^ + tp) : (3.26) 
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where A ~ o 1 Ao- 

Let us notice that the factor e lk f x in the eq.( |3.25| ) is an inherent part of of the second factor. 
In fact momentum H a (x) has a constant (x-independent) part Ilo which corresponds to a motion 
of a soliton without changing its configuration. This part is also called zero mode and has to be 
treated separately. An electron with a momentum close to the Fermi surface corresponds to a state 
with Ilo = k^. The factor e Ji ^ r S( kf ) = ± ( a relative phase between left and right movers) is more 
subtle, and not too important in one dimension. In Sec |4.3.2 we shall see how this factor develops 
into two spatial dimensions. 



3.2.3 Hydrodynamics and spin-charge separation 



In terms of the bose fields the Lagrangian (3J3) becomes 



L = - {d^f + ~(3^) 2 + A( cos(2^ T + <p) + cos(27r^ + <p)) + 

-6,^8^ + ^) (3.27) 

where we introduced an external electromagnetic field to keep track of the response functions. 
Introducing charge and spin densities 



V2 ' YS V2 

we get 



(3.28) 



L = \{d<\>cf + \(dcf) s f + Acos(V27r0 c + <p) cos(V2^0 s ) + V^e^A^d^ (3.29) 

At energy less than A the field tp follows the field —\^2tt(/) c in order to keep the argument of the 
first cos fixed, i.e., a singlet electronic current adiabatically follows the density of solitons of the 
phonon field = ^e^dyip . The dynamics of the remaining variables consists of an independent 
gapped spin density wave and a gapless charge density wave mode: 

L = L C + L S , (3.30) 



Lc = |(30 c ) 2 + V^e^A^d^c, (3.31) 

- V2iT(j) c = if, (3.32) 

L a = ^(d(f> s ) 2 + Acos(V27T(/g. (3.33) 

The charge density sector is a hydrodynamics with compressibility differs from the compressibility 
of free fermions by a factor 1/2. 

The lesson we may learn from the Frohlich's story can be summarized as follows: 

If the system has infinitely many topologically distinct degenerate vacua \Q), then a topological 
configuration (p(x, t) which transforms one vacuum \Q) into another |Q + 1) ( (p(x = ±oo, t = 0) = 
and <p(x = ±oo, t = oo) = 2tt ) is a hydrodynamic mode, regardless of whether the particle spectrum 
is gapped. Then, "the whole system, electrons and solitons, can move through the lattice without 
being disturbed" E 
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3.3 Correlation functions 



3.3.1 One particle electronic Green's function 

The bosonization approach reviewed in the previous section provides an effective way to compute 
correlation functions. Below we concentrate on the one particle Green's function 

G(t,x) = (cl(t,x)c a {0)}- 
The Green function consists of two chiral parts 

G(x,y) = e ik f^-y\^ R (x)My)) + e- tk ^ x -yH^ L (x)^ L (y)) + 

As an implication of the Orthogonality Catastrophe, the part which does not conserve momentum 
vanishes in a macroscopic system (L — > oo) 

<^(*)^(y)>-0 (3.34) 

Having in mind an extension of this method to higher dimensional topological liquids, let us 
introduce the following notations: ipk t (x) = V>l..r( x )> 0k f (x) = and G(kf,x) = Gl,r(x) 

for kj = ±kj. Sometimes we also use vj = ky/m = ±Vf. Then we can write the Green function 
as 



G(x)= £ e ik ^(4 / (x)^ k/ (0))=^G(k / ,x). (3.35) 



r «k f x/„/,t 

Bosonization decomposes the Green function into spin and charge parts 

G L (x) ~ ( e ^ (l ^) {e ^C*f W-*f W)). (3.36) 
The charge part is easy to compute, since the charge sector is just an ideal liquid 
i((^ c (k/,s) - ^(k/, 0)) 2 = ln(( V/ (t - *0) - x)k f ), 

i((^(x) + <^(0)) 2 > = ln(^( V/ (t - iO) - z)k f ). (3.37) 

The correlation function is determined by the compressibility, which differs from free fermions by 
a factor 1/2: 

(e ^(^(k / ,,)-0 c (k / ,O)) ) ^ (kf ( ^ (t _ . Q) _ x)) l/2_ (3 _ 38) 

A similar factor 

which occurs in a nondiagonal Green function ( |3.34 ) vanishes as (pr(v/(t — iO) — x) 1 / 2 . 

The spin part of the propagator can not be calculated in the same manner. However, it can be 
estimated at the mass shell \vj x 2 — t 2 |Ag 3> 1. First of all the Lorentz invariance requires that 
the Green function ( |3~35|) to be of the form (v f (t - iO) - x)~ l f{{x 2 - (v 2 f t 2 ) A 2 ,) . This gives the 
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form of the spin factor in ( ft.36| ) . It is (v/(i - iO) - x)- 1 / 2 /((x 2 - {v f t) 2 )Al). Then at t = 0, the 
spin factor in (|3.36| ) may be asymptotically replaced by the Green function of massive fermions 



where 



G°(k f ,x,t) = J ^PjP*+^G°(u,k f +p), (3.40) 

G°(u,k)= V = 17^7^2 (3-41) 

is the free fermionic Green function. As a result the mass shell asymptotic reads 

(e^<°») ~ (^) 1 / 2 ( X "^ + ^° ) 1 / 4 G° L (x,t) (3.42) 
«;/ x + Vft — i0 

so that the mass shell behavior of the Green function is 

G(kf , x, t) ~ D(k f , x, t)G°(kf , x, t) (3.43) 

where 

is the propagator of the soft modes. In a more general case, velocities of charge and spin waves 
may be different. Then 

D(x,t) ~ ( j -) 1/2 ( X "^ + '° ) 1/4 ^ Arg(kf) (3-45) 

where Arg(kj) = 0, tt is the angle of the Fermi vector. It is instructive to see the Green functions 
in the momentum representation. It is a convolution of the propagators of free massive fermions 
and a sonic mode of density modulation 

G(w,k) = £ /^(k-kz-q^-n^n)^, (3.46) 
where k and q have only two distinct directions along kj = ±/ey. Here 



1 /fr - v iq 2 + i0)\ 1/4 ,«,.„,, \ 
D(q,Q) ~ — i-5 f ^ -) e *Arg(k f ) 

is the propagator of the density wave. 

Computing this integral one gets on shell asymptotic of the Green function 

E i 2 17 v f,7 k j s = v%k = k / )2 = A ° + 'V'^ 



(kf). 

(3.47) 
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3.3.2 Electronic spectral function 



Electronic Green function ( p. 47 ) loses its poles at uj = ±y^u 2 (k — kj) 2 + A 2 ,. This happens because 
of the interaction with densities soft modes. Instead of a pole the Green function now has two branch 
cuts, starting at u) = ±yvj(k — ky) 2 + Aq. This means that an electron is no longer an elementary 
particle. It is a composit object made of the solitons of the charge and spin sectors. It is instructive 
to write the Lehmann representation for the electronic Green function 

G(4 ,„) _ P %g) + B ("f\ UE. (3.48) 



It consists of the convolution of the spectral functions of massive fermions ( p.4l|) 

A (p,E)=B (-p,E)= £ (E-v f p)5(E 2 - v)p 2 - Ag) (3.49) 

Vf=±Vf 

and the spectral function of the soft modes 

P(q,e) = —QmD(q,e) ~ e] 1/2 (e 2 - v}q 2 )- 1+1/4 ,at e > v f \q\, (3.50) 



so that 



A(k,E) = / A {k-q,E-e)P(q,e 
Jo 

B(k,E)= | o °° B (k-q,E-e)P(q,e)^. (3.51) 



dedq 
dedq 



Computing the integral over energy, we get 

-1/2 f n (k f + k-g) 

((J W 2 (k/ + k _ q) 2 + A 2 _ E) 2 _ v 2 q 



A(k, E) ~ e^ 2 — -dg, (3.52) 

y,) f In OT7T. ~ _n 2 7 A2 _ p\2 _ „2„2\l-l/4 y ' V 7 



where 



, n 1/-, Vffc 



2 .^^2 + A 2 

is the occupation number of free massive fermions. 

The integral in ( |3.52| ) goes over the branch cut, which starts from the energy of composit 
fermions E > Vfq + ^Jv 2 {kf + k — q) 2 + Aq - the energy of a massive particle plus the energy of 
soft modes. 



The eq. ( ^.52[) can be interpreted as follows. A composite particle means that there are many 
(not just one) states with a given energy and momentum. In addition to the momentum k and 
the energy E the intermediate states with one extra particle are characterized by an additional 
quantum number q — the momentum of the soft modes. The amplitude of this state \E,k,q) is 
just a free fermionic wave function *$>o(k — q), but the number of these states is 



u(E,k,q) = e] 1/2 ((^v%k f + k - q) 2 + A 2 - E) 2 - v 2 q 2 y l+1/ \ (3.53) 
so the spectral function is 

A(k, E)= J |* (k/ + k - q)\ 2 v{E, k, q)dq. (3.54) 
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The Green function is often used to compute the "occupation number" 

r°o ^E 

n(k-k f )= A(k, E) — (3.55) 

and the "density of states" 

dN(E) f „.dk , 

-L>=J A (k,E)-. (3.56) 

We have to be cautious in interpretation of these quantities. The spectral function of composite 
particles describes a tunneling process rather than a density of the spectrum. Indeed, 

A(k,E)= \(N\c(k)\\N + l,k)\ 2 5(E(N + l,k)-E (N)) (3.57) 

\JV+l,k) 

where the sum goes over all excited states of a system with M+l particles with energy E(J\T+1, k) > 
Eq(M). Here the state {j\f\ is the ground state of the system with A/" particles and E (N) is its 
energy. Were an electron be a particle (an asymptotic state), the matrix element (jV\c(k)\jV + 1, k) 
would approaches 1 at k — > kf. Then at least at the Fermi surface, the spectral function A{k,E) 
would describe the density of states A(k,E) ~ ^<5(-E(AA + l,k) — Eq(M)). In the case of the 



topological liquid, however, one particle matrix element vanishes at k — > kf (see Sec. ( 3.5.1 )). 
As a result, the spectral function reflects the matrix element rather than the spectrum. The 
spectral function appears in processes with emission and absorption of electrons, such as tunneling, 
photoemission, etc. In what follows we refer (3.55, 3.56] ) as tunneling occupation number and 
tunneling density of states. 

The tunneling occupation number is 

n(k) = / no(k — q)v{k, q)dq (3.58) 



where 

"(k,q)= I i/( J E,fc,g)^~e7 1/2 (^)- 1 / 2 atg»A (3.59) 

where the last factor in the integrand is the number on states with a given momentum. 

The tunneling occupation number has a broad character of a scale kf around the Fermi points. 
At the region kf ^> \k — kf \ S> Ao, where the gap is not important we get a familiar result for its 
singular part 

n(k) - n{k f ) canst |Jfe/Jfe/ - l| 1/2 sgn(/c - kf). (3.60) 

It is a smooth function at k = kf and crosses kf linearly: n(k) — n(kf) ~ — const (| A; — kf) at 
\k — kf \ <C Ao- Contrary to the free massive particle case, the tunneling occupation number does 
not vanish at the Fermi surface. This does not mean that states with k > kf are occupied in the 
ground state. They become occupied in the process of embedding an additional electron into the 
system. 

The tunneling density of states behaves similarly 

d ^=lAik,E) d ±~e-^( E dN °V (3.61) 
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where dNo/de = e/^J e 2 — Aq is the density of states of free fermions. The density of states shows 
an asymmetric broad (~ £/) peak. It decays from the peak toward the threshold E = Aq as 

dN }^ ~ const |£/ e/ | 1/2 , at e f > E - A > A (3.62) 
dE 

In contrast to free massive particles, the tunneling density of states is not singular at E = Aq. It 
remains smooth at the threshold and approaches it linearly (~ E — Aq). Among the interesting 
features of the spectral function, one is worth mentioning: the spectral function is determined by 
two scales, Ao and e/, rather than just Ao- The second scale is the signature of the orthogonality 
catastrophe and of the composite nature of the electron. 



3.3.3 Axial current anomaly 

One can reproduce the result of the last section by performing an anomalous transformation 

i>L* = e^^XLa, jilt* = e-&*XL*, (3-63) 

where pl and p>R are holomorphic and antiholomorphic parts of the field p, so that (pt±v td^(pL,R = 
and pl + Pr = 2<P- Then formally the field (p may be absorbed by the chiral gauge transformation 
■^lr ~^ -A^lr + ^xPl,r so that it disappears from the Lagrangian ( |3.9| ). This is, however, not true 
due to the axial current anomaly. To correct it one must add to the Lagrangian the hydrodynamics 
of zero modes, namely the term ^2(|^) 2 , and connect it with charge density ( |3.12| , |3.32| ). One 
can see this immediately by comparing a bosonized version of the Lagrangian before and after the 
anomalous transformation. The bosonization treats an anomalous axial symmetry correctly. 

The operator Xa carries no charge [Xa^c] = and is a soliton operator of the spin sector. 
Due to the constraint ( |3.32|) , its bosonized version is xt ~ oT 1 ! 2 : e^^ 3 ^ 1 '^ : and xi ~ aT 1 ! 2 : 
e v?^ 8 ^!'^ : . In its turn the vertex operator 

V c (k f ,x) = e^AxgO^Mk,,*): (3 . 64) 

annihilates a soliton in the charge sector. 

In terms of %a particles, the Green function has the form 

G(x)=J2(xi f Jx)V c 7 1 (k f ,x)V c (k f ,0) Xkf A0)), 
so the mass shell asymptotic of the Green function (p. 46) becomes transparent. At equal times, 



where there is no complications with the Lorentz factor (|3.42|) one may treat charge and spin sectors 
independently 

G{x) = J2(xl } »Xk /lCT (O))£0z) (3.65) 
k/ 

where 

D(x) = (y~ 1 (k / ,x)y c (k / ,0)) (3.66) 
is given by fl3.44|) . This reproduces the result of eq. ( 0.43[) . 
Under the constraint (3.32) the vertex operator becomes 



V c (k f ,x)=e^ Al ^ )M ^ xy -. (3.67) 



It carries no spin but annihilates a unit charge [V c (x), p(y)] = V c (x)5(x — y). 
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3.3.4 Pair correlation function and a long range order 

Pair correlation function. Among many possible two-particle correlation functions, let us con- 
sider the pair correlation function 

K(ri,r 2 ,r 3 ,r 4 ) = (cj(ri)cj(r 2 )q-(r 3 )q(r 4 )) (3.68) 



It illustrates a difference between the tunneling amplitude (|2.23| ), the order parameter and the long 
range order. 

The function consists of harmonics e lk f{ T ^ r "i- r ^ T i) . Let us pick up one harmonic, say, e ik f( ri ~ r2 ~ T3+r4 \ 
which corresponds to a correlation between pairs located at r\,r 2 and r 3 ,r 4 : 

K(n,r 2 ,r 3 ,n) = e^^-^-^^^Cri)^^)^^)^^)}. (3.69) 
Bosonization gives 

K(r 1 ,r 2 ,r 3 ,r 4 ) = a - 2 e % (ri_r2 ~ r3+r4) ^(^R^-^lA^)) ^^{^rA^-^lA^))^ 

^ e ^(<t>R, s (ri)+<t>L, s (r2)) e -^{<t>R,s(r 3 )+4> LtS (r 4 )) , ^ ^ 

Let us assume that the size of a pair is less than the distance between pairs \r% — r 2 \ ~ \r 3 — r 4 | <C 
R = \r\ + r 2 — r 3 — r 4 \ and let us consider the dependence of the correlation function on R. We 
have 

K(R) = (^(R)^ Ll (R)^m Ll (0)) ~ a- 2 (e^r-W^)( e — MMR)-Mo)) } (3 . 71) 
where 2-7T0 = 4>r + (f>L, 2ir = d x ((pR — <Pl)- The spin sector factor approaches to a constant at 

R — > OO 

/ e -iirV2(^ a (R)-(j> a (0))\ ^ -y (3.72) 



whereas the charge sector factor can be computed by means of eqs. ( 3.31 , 3.32 ) 



K{R) ~ ^ (3.73) 

The correlation between pairs decays, but not as fast as a similar correlation for free particles 
(~ R~ 2 ). A decay of the pair correlation function must not be misleadingly considered to be a 
signature of the absence of superfluidity in one dimension. As we discussed above, superfluidity 
indeed does not exist in one dimension, but for a different reason. The reason is a pinning by 
impurities, rather than a decay of the correlation function ( |3.71 ). In the next sections we show 



that in spite of the decay of K(R) not only the superfluid density, but also the matrix element 
(ct(ri)ct (r 2 )) does not vanish. 

A long range order. It is easy to construct an object which carries the electric charge 2e and 
exhibits a long range order. For instance, the operator &(ri,r 2 ) = ct(ri)V c (ri)V c (r 2 )(y,(r 2 ) shows 
the long range order: 

K{R) = (*t( ri ,r 2 )*(r 3 ,r4)) ~ 

( e iVl(<pRA r i)+<pLA r 2)) e -iV2(<l> R Ar3)+<i>L,s(r4))} ^ I (3.74) 

at R — > oo. This long range order indicates the Peierls instability at 2k t. 
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3.4 Vertex operators and local singlet state 



3.4.1 Vertex operators and operator algebra 

In this section we rederive the results for the Green function (obtained by the bosonization in the 
previous section) by means of the Hamiltonian approach. This approach, although less automatic 
than bosonization, has certain advantages: it clarifies the physical picture and can be lifted to 
higher dimensions. 

Let us introduce a vertex operator of the spin sector 

V s (k f ,x)=e^ :Mkf ' x) \ (3.75) 
This operator is charge neutral [</> c > Vs] = but creates spin 1/2 state: 

[S 3 (x), V s (y)} = -\v s {x)5{x - y) (3.76) 

[S 3 (x), (VM)- 1 } = \{V s (x))-H{x - y) (3.77) 
where S 3 (x) = -k=d x (b s is a spin density. The charge vertex operator, then, as a composition of 

V ^7T 

the electron and an antikink of the spin density. The latter removes the spin from the electron 
operator. 

The vertex operator obeys an operator algebra. 

h/V )k », V2^(k /9 0)] = °H-j^) (3-78) 

and therefore 

1>k f Ax)V.(y) = ^(y)^k /lCT (^)|^r /2 e^ Ar S (k ^ (3.79) 

where Arg(ky) = 0, tt is the angle of the Fermi vector and a = ±. Being written in this form the 
equations ( |3.79j ) have a straightforward generalization to higher dimensions. 



3.4.2 Local singlet state 



Let us consider a system with periodic boundary conditions and even number N of particles, such 
that the ground state is a singlet \N). Now let us attempt to embed an extra spin up electron. 
According to ( |3.12 ) a state with one extra particle requires a change of boundary condition ( one 
extra particle requires a half of a kink <f>{x = oo) — cj)(x = — oo) = tt) from periodic to antiperiodic. 
Also, the spin of a state with odd number of particles and antiperiodic boundary conditions will 
be 1/2. Technically it will be convenient to consider the one particle tunneling process by keeping 
the boundary conditions periodic. One has to add two particles (compatible with the boundary 
conditions) but separate them in order to consider the particles independently. The new state 
remains a singlet. Moreover, the spin localized in the vicinity of each electron bound to a soliton 
is also zero. Following Laughlin we refer it as a local singlet. 

How does one create the local singlet state \M + 1)? One can try to act with the operator 
V>j(ky,x). The state ipt(kf, x)\J\f), however, has spin 1/2. To create a local singlet state we must 
remove the spin from the electron by creating a proper 1/2 vortex of spin density acting by the 
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vertex operator V s . To do this, let us first write the electronic operator Xa in terms of modes of 
free massive Dirac particles 

Xa (x) = J (^ ) (k,x)a (T (k) + *W(k,x)ot (k))dfc (3.80) 

where 

<4 +) (k,r) = u k _ k / kr + v k _ kf e^- 2k f>, (3.81) 

cfrjf^r) = v k „ kf ^ - u k . kf e^~ 2 ^ r . (3.82) 



are Dirac wave functions of positive (negative) energy E p = Jp 2 + Aq and u p and v p are 



1(1 + ^) (3-83) 

5(1 (3-84) 

In these terms we can write the operator which annihilates (creates) a singlet state with one 
extra particle with momentum p relative to the Fermi surface is a composition of an electron and 
the vertex operator 

a(k f ,x) = a T (k /5 x) : V s (k f ,x) := ai(k f ,x) : (K(k/, a;)) -1 : 

~(%) 1/2 :V s (k f ,x)a ] (k f ,x): (3.85) 
The vertex operator in (3.85) removes the spin from the electron by creating a soliton and 



binding the soliton and electron. One can reinterpret the eq. ( |3.85| ): an electronic operator is a 
composition of two operators which create half of a soliton in the spin and charge sectors. The 
vertex operator annihilates the spin-half-soliton. Therefore the operator a creates one half of the 
charge soliton and inserts a particle. This gives the true ground state with one extra particle: 

\M+ l,k f + q> ~ J a\kf,r)e- l( ^f + ^ r dv\N) (3.86) 
It carries a unit charge but no spin. 

3.5 Matrix elements 

The operator algebra and the constructed local singlet state allow one to compute matrix elements. 
We start with the one-particle matrix element. 

3.5.1 One particle amplitude and Orthogonality Catastrophe 



Using eqs,( [3.85" , 3.86 ) we can write the one particle matrix element as a correlation function in the 



ground state with J\f particles. 

(A/'|c T (x)|A/'+l,k / + q) ~ 

J (|) 1/2 (AA|^ k/ (x)y;(k / ,y)aj C/iT (y)|AA)e ik /( ;c -^-^dy ~ (3.87) 
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Here we used the operator algebra of the Sec. 3.4.1. The matrix element of the spin vertex operator 
on the ground state is (J\f\V s (0)\N) = 1, and we have 

(Af\c^k)\AT + l,k) ~ e^SO*/) f $< +) (k/, k - q)v{q)dq (3.88) 

where u(k,q) is the number of states with a given momentum ([3.59 ) . 

Similarly, the wave function of the one particle eigenstate, i.e., the one particle matrix element 
with a given energy and momentum is 

(Af\ Cl (k)\Af + l,k,E) ^fArgCk/) f $W(k / ,fc - q)u(q,k,E)dq (3.89) 

where z^(q, fc, E) is given by the eq.(p3|). 

The overlap between the wave functions of the ground states of system with M and M + 1 
particles suppressed by the factor \k/kf — l] 1 / 2 at \k — kf \ 3> Ao- At \k — kf \ <C Ao the matrix 

1 /2 

element, is nonzero, but small ~ (Ao/ef) 1 '. This ph enomenon is known as an orthogonality 
catastrophe. A physical reason for this is that a spectral flow is always even degenerated. A 
state with one extra particle fractionally (1/2) occupies the zero mode created by a soliton only 
fractionally 1/2. Any fractionally occupied state is orthogonal to fully occupied states. 

Similar arguments may be applied to any odd particle amplitudes. They vanish as a power of 
k — kf the size of the system 1/L or a gap A, whatever is smaller. 

Matrix elements between states with an even number of particles are different. Low energy 
states with even number of particles are nondegenerate (completely filled) — an overlap between 
nondegenerate states is nonzero. 



It is instructive to rewrite the matrix element ( p.88|) in different forms. It can be written as 
an expectation value over the ground sate with M + 1 particles of the vertex operator of the spin 
sector: 

(N\ci (k)\N + l,k) -(a/L) 1 / 2 J dxe ikx ^\k f ,x)(Ar+l,k\V s \k f ,x)\M+l,k} (3.90) 

Using the operator algebra, we have 

(AA+l,A;|F;(k / ,x)|AA + l,fc)~ e fc (^) 1 /2(AA|y;(k / ,x)|AA) (3.91) 

The latter matrix element is 1. This shows that the vertex operators have different matrix elements 
on states with odd and even number of particles. 

Another instructive form is found by virtue of the vertex operator of the charge sector. Accord- 
ing to flgnnyrgap 

c a {x) =V- 1 (k f ,x)xAk f ,x) (3.92) 

so, we have 

(Af\c a (k)\M + l,k) ~ 

J dxe lkx (Af\V c (k f ,x)V c \k f ,0)\Af)(Af\xk f A^+^) (3-93) 

where the matrix element of operator Xa is given by eq. ( |3.81| ), and the correlator of charge vertex 
operators is (3.44). 
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3.5.2 Two particle matrix element 

Let us use the same technique to compute a two particle matrix element: 

A(x-y) = e (Ta/ (N\c fT (x)c a ,(y)\Ar + 2) (3.94) 

As we shall see the result will be different from the one particle case. There are three related 
differences. First of all one can embed two particles with a zero total momentum by putting them 
on the opposite sides of the Fermi surface. Secondly, a two particle state can be a singlet even 
without the help of solitons. Finally, by adding two particles we create and fill completely the zero 
mode. As a result the orthogonality catastrophe would not show up for the two particle case. 

The two particle matrix element is 

A(x-y) = ^(AA|c T (x)c i (y)at(k / ,x)at(-k / ,y)|AA)e ik /^ (3.95) 



k 



/ 



or, using ( |3.85j ) 

A(x - y) ~ E k/=±fc/ (AA| (^~ y H L ^x)i> R ,i{y) + e-^-^ jT (x)V>L,|(y)) 

V s {k f ,x)V- l (-k f ,y)a){k f ,x)a){-k f ,y)\M) (3.96) 
Proceeding the same way as in the previous section by using the operator algebra ( 3.7S| ), we find 



{kfx) 1 / 2 

where 



. . . sin(kfx) . , . ,„ 
A ( x ) ~ n \i/9 ^o(x (3-97) 



A (x) = v f f sign(A; - k f )$ { +) (k,x)<I> ( +) *{k,x)dk = v f f e ipx =dp (3.98) 

■i ■' lv 2 f p 2 + Ag 



is the matrix element for free massive particles. Notice that the latter has the form of the BCS 
wave function. In the one dimension it is the modified Bessel function 

A (x) = A K (A x). (3.99) 

This formula gives the universal part of the matrix element as a function of Vf(k — kj)/Ao at 
\k — kf \ <C kf : 



^)-A( fc/ ) „ ( } { const + (^)V2 l n (^) Vfp » A 



The rest of the function, including its value at the Fermi surface Ao(&/) ~ (Aoej) 1 / 2 is not universal 
and is determined by processes with large momentum transfer. A qualitative graph of A(k) has 
a broad non-universal asymmetric peak deep inside the Fermi-surface and approaches the Fermi 



surface according to ( 3.100 ). At \k — kf\ ~ Ao, it changes behaviour and crosses the Fermi surface 
linearly . As we can see A (A;) is drastically different from Ao(k). 

Let us comment that there is no contradiction in the fact that the correlation function K(R) 



(^(R^l^R^R^O^niO)) vanishes (g7| gjj ), while the matrix element A(R) = {^(R^l^R)} 



27 



does not. The reason is that the "clusterization theorem" which connects the long range order and 
the matrix element is valid only for a true condensed system where all states with a nonzero mo- 
mentum are separated from the ground state by a gap or for a system where the matrix element 
between an electron and the soft modes is small, like in BCS. Neither is the case in our system. In 
the Lehmann expansion of 

K ( R ) = E E elPR \ Wtr (°) W + *,frP) 1 2 
p </> 

the are many intermediate states with the same given arbitrary momentum P. They are character- 
ized by the profile of solitons <f>(x), i.e. by configuration of soft modes, and have a singular density, 



similar to ( 3.53 ). A sum of oscillating factors decays with R. 



3.6 Concluding remarks 

As we have seen the Frohlich-Peierls models provide all necessary spectral features of a supercon- 
ductor: 

• the electronic liquid is compressible - the superfluid density is a smooth function of chemical 
potential; 

• one particle spectrum is gapped; 

• a pair wave function is localized and normalizable; 

• there is a long range order, although the order parameter differs from (ctcj). 

However, Frohlich's electronic state is not a superconductor, but rather an ideal conductor, because 
of the strength of pinning in one dimension - an arbitrary small concentration of disorder generates 
a resistivity. 

Topological liquids in higher dimensions discussed in the next section are free from this flaw. 
Pinning plays no role as long the system is a liquid. 



4 Topological superconductivity- Two Dimensions 

In this section we discuss the simplest model of the topological fluid in two spatial dimensions. 
Its physics is very close to the one dimensional example of the previous section. To stress the 
similarity, we try to follow the line of Sec. 3 as much as possible. The two dimensional model is not 
as well developed as the one dimensional model. In particular the identification of fast and slow 
modes and the passage to a continuum field theory from a microscopic lattice model are not as 
straightforward as in the case of the Peierls model. We therefore start directly from the continuum 
model. For the account of attempts to derive this model from the doped Mott insulator, see |2j. 

4.1 Topological liquids and topological instability 

Let us consider an electronic liquid where the interaction between electrons is mediated by an 
electrically neutral bosonic field, that can form a point-like spatial topological configuration (a 
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soliton). Let us suppose that in a sector with zero topological charge the electronic spectrum has 
a gap Ao- Assume now that the system exhibit a spectral flow. This means that in the presence 
of a static soliton the electronic spectrum differs from an unperturbed one by an additional state 
just at the top of the valence band or within the gap — a so-called zero mode or a midgap state 



[19]. If the zero mode is separated from the spectrum, its wave function is localized around the 
core of the topological defect. In case when the level is attached to a band, the wave function 
decays as a power law away from the center of the soliton. A general argument |2(J suggests that 
the midgap state always has an even degeneracy. This degeneracy eventually leads to a proper flux 
quantization and below it is assumed to be twofold. 

Now let us add an even number of extra electrons with the concentration 5. They may occupy 
a new state at the Fermi level of the conduction band. It costs the energy of the gap plus the 
Fermi energy Ao + [i per particle, where [i is a chemical potential. Alternatively, the system may 
create a topological configuration and a number of zero modes in order to accommodate all extra 
particles. The energy of this state is the soliton mass plus exponentially small corrections due to 
the interactions between zero modes. If the latter energy is less than p, + Ao then every two extra 
electrons added to the system create a soliton and then completely fill a zero mode, rather than 
occupy the Fermi level of the state with zero topological charge. As a result the total number of 
solitons in the ground state is equal to half of the total number of electrons in the system. 

Formally it means that, contrary to the Landau Fermi-liquid picture, the expansion of the 
energy in a small smooth variation of chemical potential 6fj,(r) with a fixed topological density of 
solitons F(r)/2ir has a linear term in <fyi(r): 

5E(jj) = - J 6n(r)p(r) dr + J 5n{r)K{n, r - drdv' + 0{5p 2 ). (4.1) 

with a non vanishing zero harmonic of the kernel J K(p,r)dr ^ 0. The linear term in chemical 
potential is known as Chern-Simons term. 

The minimum amount of energy is achieved if the variation of density follows the variation of 
the topological charge 

p(r) ~Ps = j K(p,r- r')^p-dr'. (4.2) 

On doping, electrons create and occupy zero mode states to minimize their energy, thus giving a 
non zero value to the topological charge (compare with (|3.4|) ): 



5={J K(^v)dr) J ^dr. (4.3) 

This is a topological instability. 

This is already sufficient to conclude about superconductivity — the chemical potential always 
lies in a gap. The arguments are the same as in one dimensional case and can be borrowed directly 



from Frohlich's paper [14]. The position of the topological excitation is not fixed relative to the 
crystal lattice. Therefore, a pair of electrons bound to a topological excitation can easily slide 
through the system (and carry an electric current). It slides without attenuation, since the state 
is completely filled and is separated by the gap from unoccupied electronic states. As a result, the 
low energy physics of density fluctuations is described by the hydrodynamics of a liquid of zero 
modes (2/7). Indeed, due to a gap, the system has a finite rigidity to a local change of the density 
of topological charge: 

5E = £^ J F 2 {r)dr 
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. This together with ( [4.2j ) yields the compressibility (|2.7| ) and the hydrodynamics. In dimensions 
higher than one the hydrodynamics already implies the Meissner effect (2^,2^) and superconduc- 
tivity, g 

A standard example of the theory which exhibits both spectral flow and a topological instability 
and therefore superconductivity is the Dirac Hamiltonian with a fixed chemical potential [] 

H = £.=1,2 " + A)t/> ct + A /3 + MO))^, (4.4) 

Here a = (a x ,a y ) and (3 are 2x2 Dirac matrices: {a x ,a y } = 0, /3 = —ia x a y and the gauge field 
A mediates the interaction between electrons (it is not an electromagnetic field). For a relation 
between this model and doped Mott insulator see Q. 

This Hamiltonian has / F(r)dr zero modes (the flux is directed up). Wave functions of zero 
modes are: 

$(r) = e~ l f Mr'ydr'-pf A(r>)xdr' M -^ 

0^ o (z) = -Q Q (z), 
aV$ = 0, 

where <3?o is any polynomial of degree ^— J F(r)dr — 1 [24]. If the soliton has a unit topological 
charge <l?o is a constant. 

This single fact implies that the energy at = has a linear term (t4,l|) in chemical potential 



j K{jx = 0,r)dr = 2 



. There are other models of topological superconductivity in any dimensions (Sec.l), which we do 
not discuss here. 

Let us comment on the relation between the topological mechanism of superconductivity and 
superconductivity in a system of anyons [16, [l7|]. The models become very similar after projection 
onto the low energy sector. Then the relation fl4.2| ) can be imposed as a constraint rather than as 
a result of minimization of the energy. The projection can be done by introducing a Lagrangian 
multiplier Aq for the relation ( |4.2| ) and commutation relations 

[A x (r), A y (r')] = 2TTiK-\v - r'), (4.5) 

or by adding the Chern-Simons term with the kernel K to the Lagrangian. Also at large gap and a 
small concentration S, one may replace the Dirac Hamiltonian with the Pauli Hamiltonian, so that 

H = ^-(-iV + A) 2 + A (4.6) 



In anyon model the kernel K(r) is replaced by 25 (r) (even at /j, ^ 0), so that the relation 
between topological charge and the density ( f4.2j ) becomes local p(r) = 2-^f~- This simplification 

4 This mechanism is very different from a condensation of pairs bound to a polaron. fn contrast, in the topological 
mechanism the electric charge of a pair is only partially localized at the core of the topological soliton. Although the 
number of zero modes is equal to the topological charge of the soliton, a part of the charge is smoothly distributed 
throughout the rest of the system. 

5 The same model, written with a fixed number of particles looks different. It requires an additional commutation 
[Ai(r), Aj(0)] = 2ni€ijK~ 1 (r) or a proper Chern-Simons term in the Lagrangian formulation. 
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results in a generation of transversal electric currents or an "internal" magnetic field by light or by 
inhomogeneous electric charge. For the same reason, the Meissner effect in the anyon model per se 
exists only at zero frequency, zero momentum, zero temperature, infinitesimal magnetic field etc. 
An origin of these unphysical consequences is the topological constraint ( |4.2| ) is taken locally and 
instantaneously. In topological superconductors the kernel K(p, r) is determined self-consistently. 



4.2 Hydrodynamics 

4.2.1 Compressible charge liquid 

The hydrodynamics of a superfluid ( 2.161 ) can easily be obtained from the model (4.4). Let us 
see how the energy ( |4.lD of a spin singlet state changes under smooth variations of "electric" 
E = dtA — and "magnetic" F = V x A fields. To keep track of spin variations we add 
an external field A 3 to the Hamiltonian ( |4.4[ ) iV — A — > iV — A — .A 3 03. To obtain a linear 
hydrodynamics it is sufficient to compute radiative corrections in the Gaussian approximation. Let 
us now fix the chemical potential and set the total spin to be zero. The result consists of two 
separate parts — the spin sector and the charge sector: 

5E = 6E C + SE S , (4.7) 

In the Coulomb gauge VA = the density of energy is 

SE C = ^(vftf + F^+S^K^-Sp), (4.8) 

SE S = ?±( v f(d t A 3 ) 2 + (VxA 3 ) 2 )+^A 3 VxA 3 , (4.9) 
where polarization operators 

k 2 U± {u,k) = (u, k)j ± (-uj, -k) - j\\ (u, k)j\\ (-u, -k)) = 

Kl^ ~ 5 ll )jTra l G (p)a,G (k-p) J ^, (4.10) 
u;K(uj,k) = i(j l] (uj,k)i ± (-u J ,-k) = ie ij jTta l Go(p)a J Go(k-p) J ^ (4.11) 

are current-current correlators (j = jj_ + jn) of free Dirac massive fermions and Gq(p) = (pq — ap — 
/3Aq) _1 is the Green function of massive Dirac particles. 



Minimization over 5p gives the relations (4.2) for the charge sector. Substituting this relation 



into the rest of (4.8) we get 



,2vr N 



6E c = QfU ± ((8p) 2 +vfi 2 ) (4.12) 



This gives the hydrodynamics of the ideal liquid ( |2.7| ) provided that the propagators (2tt) 2 K 2 (oj, k)H±(u 
are nonzero and analytical at k, lo — > 0. 

The latter is indeed true. The off-diagonal propagator K(uj, k) fixes the number of zero modes 
of Dirac operator per unit flux, so K(0) = 2, whereas H±(uj,k) describes vacuum polarization of 
gapped particles. Some caution and proper regularization is necessary, since the integrals diverge 
at large p, but otherwise the calculations are standard 

Il ± ~£ o (l + O(k 2 /A 2 )) (4.13) 
K = 2 + 0(k 2 /A 2 ). (4.14) 
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4.2.2 Incompressible chiral spin liquid 



In the spin sector (£0j) the story is different. The total spin is kept to be zero. Therefore, the Chern- 
Simons term (the last term in ( |4.9| )) remains in the spin sector and results in the hydrodynamics 
of an incompressible spin liquid. Writing <9fU s = — 5C/5A 3 , we obtain the hydrodynamics of a 
topological spin liquid: 

C s = 2vru s x d t u s . (4.15) 



By combining (4.8, |4.9| ) we obtain the hydrodynamics of the topological superconductor: 

The hydrodynamics consists of two independent fluids C = C c + C s : a compressible charged 
liquid 



k), vfi (— u, — k)) = k 



-1 



v 



k2' 



(u ± (w,k),u ± (-w,-k)) = kT 1 -^, (4.16) 



and an incompressible topological (chiral) spin liquid: 

(ul(u,k),u^(-uj,-k)) 



Attlo ' 



where u", u are the longitudinal and transversal parts of the displacement u«(k) = -£u" H — ^-u . 
Eqs. ( [4.16 ) imply compressibility of the charge fluid and the Meissner effect. The hydrodynamics 



d4jj) of the spin sector is equivalent to the hydrodynamics of the FQHE fluid (see e.g., pi] 



One of the direct consequences of the chiral nature of the spin liquid is that the spin liquid 
generates spin edge current. Similar to the FQHE the spin excitations are suppressed in the bulk 
but develop a spontaneous spin edge current with the level k = 2 current algebra. Let us stress that 
the spin edge current is the only hydro dynamical manifestation of spontaneous parity breaking. 
Contrary to a number of claims scattered through the literature, the spontaneous parity breaking is 
invisible in the charge sector even for the systems with an odd number of layers, but is found in the 
spin sector in a system with a boundary. A similar phenomenon takes place in the one dimensional 
Peierls model with open ends. In this case there are free spin excitations on the ends, although the 
bulk spin excitations are gapped. 



4.2.3 London penetration depth 



Equation (4.14) gives the scale of compressibility and the London penetration length. From 



h 2 v 2 f /A . This result is very different from the conventional BCS the- 



Q4.12jj4.14D we have f 

ory. In the BCS the compressibility and the penetration length are determined by the plasma 
frequency ( a classical object, independent of fi and the gap in the spectrum). It is given by the 
atomic parameters and the density of conducting electrons A° = c/uj p . In contrast the compress- 
ibility of the topological liquid described by the model (44) is determined by the gap and has a 
quantum origin. The reason for this difference is the same as for the orthogonality catastrophe - 
conducting electrons are not quasiparticles. Only a small fraction of the conducting electrons are 
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involved in the supercurrent p s /po ~ (An/e/) l , where po ~ (kf/h) d is the Fermi density (the 
volume of the Fermi sphere) and d is a spatial dimension^. 

Bearing in mind that the gap also determines the correlation length £ ~ hvf/Ao, we may 
express the penetration depth through the correlation length, atomic parameters and an electronic 
wave length £ = 2irh/kf. We write them for spatial dimensions d = 2,3. 



and the Ginzburg-Landau parameter is 

Al - A °4)^, (4.18) 



s « 

where A° = c/cjp. The Ginzburg-Landau parameter becomes universal (independent of a gap) at 
d = 3: 

^ = ^ = (c/v f )^ 2 (47Te 2 /hcr 1/2 ~ 0.1 (c/r;/) 1 / 2 ~ 10 2 . (4.19) 
Topological superconductors are of the London type. 



4.3 Electron as a composite object 
4.3.1 Electronic operator 

The most difficult part of the theory is to determine a relation between the true electron operator 
c a and the "spinned" fermion tp a of (4.4). 



There are several reasons why ip differs from a physical electronic operator. First of all elec- 
tronic states are gauge invariant, whereas ip is not. Moreover, without a gauge field (i.e. without 
interaction) tp a being a Dirac spinor has 1/2 - orbital momentum and is double valued. Electronic 
states are single valued and must have an integer orbital momentum. Another problem is that an 
electronic excitation carries a typical momentum of the order of kf, while typical momenta of Dirac 
particles are close to zero. These difficulties do not occur while studying the hydrodynamics, but 
arise in matrix elements. 

Below we conjecture the form of the electronic operator based on plausible physical arguments. 
In fact this relation can be derived (under certain assumptions) from a microscopical model of the 
doped Mott insulator (see e.g., Q). We do not discuss it here . 

The requirements for the electronic states are 

(i) the electron is gauge invariant ( in respect to mediating gauge field), i.e. remains unchanged 
under a non-singular gradient transformation A — ► A + VA; 

(ii) In a sector with completely filled zero modes, i.e., when the flux and the number of particles 



obey the topological constraint (4.2), a charged singlet excitation is a spatial scalar, i.e. its wave 
function has a zero orbital moment I = 0; 

(iii) Since an electronic liquid is compressible, the most essential electronic modes have momenta 
k ~ (2vr/5) 1 /2 = k f . 



6 This result is specific to the model (4.4). Other models of topological superconductors, describing commensurate 
doped insulator, like the model of jjj, give a superfluid density varying in the range from 5 = po (a doping) to 
pa ~ (A /e/) d_1 . 
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We find the electronic operator in three steps. The first step is the use of rotational invariance. 
Let us consider a wave function of the free Dirac field in two spatial dimensions. In the basis, where 
a-matrices are a x = 03, a y = — 02, /3 = a%, the solution to the Dirac equation with momentum p 



is 



; ipr e -i/3Arg(p) ^ ^ , E = +E p = yj p 2 + Ag 



(4.20) 



for the positive energy, and 



e ipr e -i/3Arg(p) \ E = _ Ep = _^ p 2 + A 2 (421) 

for the negative energy, where again Arg(p) is an angle of the momentum p, relative to the x-axis 
and 



3(1 - jg) (4-23) 

We choose the basis such that u p and v p to be the BCS wave functions and agree with one dimen- 
sional case (3.83). 



The spinor carries an angular momentum I = 1/2. We unwind the Dirac field by a chiral 
rotation 

Vv(p) - e^ Arg(p) ^(p). (4.24) 

This singular transformation has a clear physical sense — it projects the spinor wave function onto 
a direction of the momentum p. Indeed, the chiral transformation (|4.24j ) in two spatial dimensions, 
where f3 = —i/2 [a x , a y ] is equivalent to a spatial rotation of the momentum p by the angle Arg(p) 
which aligns the momentum p along the x-axis of the coordinate system. Now, without topological 
gauge fluctuations, the transformed operator is a spatial scalar. 

The next step is to boost the fermion to the Fermi surface. In the chosen basis the upper and the 
lower components of the Dirac field ip a = (ifi£\ w J correspond to the states propagating forward 
and backward along the vector p. To construct an electronic operator we shift the momentum of the 

upper component by the Fermi vector kf = kt — directed along the momentum p: p — > k = kf + p 

v 

and the momentum of the lower component by — kf : p — > k — 2kf = — kf + p 

( c CT (k f + p) \ f/3Ar g (k) , , s _ ( cos + ism ^|^Va 2) \ u 

l^-kf + p) )~ 62 MP) -[ism^^ + cos^^ )■ (425) 

Here we used c a to indicate that the gauge field has not been taken into account yet. 



Formally, the chiral rotation (4.24]) may be also understood in the following way. The wave 



function of the Dirac field depends explicitly on the choice of a-matrices, i.e. on the choice of 
holomorphic coordinates in a plane. The chiral transformation aligns holomorphic coordinates 
relative to each point of the "Fermi surface", i.e. sets up the momentum dependent a-matrices 

ttkf =e |/3Arg(k f ) ae -^Arg(k f )_ 

In terms of c the free Dirac Hamiltonian describes isotropic backward scattering 

H = J {64(k)c CT (k) + A (4(k)^(k - 2k f ) + h.c.) } dk, (4.26) 
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where = Vf(k — kf). This clarifies the origin of the model. It arises as the result of 2k f instability 
of some correlated electronic system. 

Let us now consider solution of the Dirac equation in the presence of the 1/2-vortex located at 

ro: 

M = ~%7^% (4.27) 
2 [r — r ) z 

The wave function in ip representation is 

/ 1 \ 1 



V -l ; (z- Zo y/^ 

where z = x + iy. This soliton describes a single particle state at the chemical potential. It is 
double valued and can not be treated as an electronic wave function. To construct the electronic 
wave function we must first write it in the momentum representation: 

where 

z p = p • r + ip x r (4.29) 

is the holomorphic coordinate relative to the vector p, and then make a chiral rotation ( |4.24| , |4.25[ ). 
This gives 

/ i(k-k^)z 
/2 dz. (4.30) 

Now the wave function is single valued, translational invariant but acquires a spatial orbital 
moment I = 1. We conjecture that it is the wave function of electronic zero mode. 

4.3.2 Bosonization in two dimensions 

The results of the previous section can be written in terms of bosonization. Let us introduce a 
chiral bosonic field <I>(k, r) similar to (|3.21|): 



$ ff (k, r), <& CT (k, 0)- : * ff (k, r), $ CT (k, 0) := In (£), (4.31) 

where z^ is given by ( |4.29| ). 

This field is consistent with Chern-Simons commutation relations and the topological constraint 

(EDQ 

d Zk $ a (k,r) = -2mu a (4.32) 

[u CT (r) x U(J (0)] = »5(r), (4.33) 

where 

u a = k • u CT — ik x u CT (4.34) 



7 Let us notice that commutation relations between vector potentials and displacements depend on whether we 
consider the theory at fixed number of particles or fixed chemical potential. In the lat ter ca se displacements commute. 
If the number of particles is fixed they do not, and the commutator is given by eq. (fh33). 
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As in one dimensional case the operator 

l/ CT (k,r) =: e* CT ( k ' r ) : (4.35) 



correctly reproduces a Dirac fermion in the sector of zero modes (|4.28| ). 

Assembling all pieces we obtain the bosonized version of physical electron in the sector of zero 
modes 

C a (r) ~ J e -*kr e -iArg(k) ^ ( k) ( k ) dk ( 4 . 36 ) 

where ^(k) is the operator which creates a zero mode with momentum k. 
Eqs. (gH-H) 

are the bosonization formulas in two spatial dimension. They express the elec- 



tronic operator through its displacement. These formulas are to be compared with the bosonization 
formula ( 3.25| ) in one dimension. 



The vertex operator V^(k, r) has a simple meaning. As is in one dimensions [p3f | it creates a 
flux quantum and a zero mode in the state with the spin 1/2. It is gauge invariant and on the 
subspace of zero modes it obeys two relations 

V r " 1 a(-iV + A) V a = -iaV, (4.37) 
[F(r'), V a (r)\ = 27rK(r)5(r - r'). (4.38) 

4.3.3 Vertex operators 

As in one dimensional case we will need an operator algebra for the vertex operator V a and two 
additional vertex operators of the charge and spin sectors. Let 

d> c =i($ T + $J, (4.39) 

$ s = I(* T -* ; ) (4.40) 

and 

V c =: e* c : (4.41) 

V s = : e* s : . (4.42) 

The vertex operator of the spin sector V s (z) creates a soliton of the spin displacement u s and 
removes spin down from the site z: 

[V s (z),Vu s (z')} = V s (z)5(z - z'). (4.43) 

In contrast, the vertex operator of the charge channel V c {z) creates a flux of the gauge field (the 
spin chirality), but commutes with displacements 

[F(z'),V c (z)} = 2ttV c (z)5(z' - z), 

[u(z'),V c (z)} = 0. (4.44) 

Due to the commutation relations (|4.33j ) vertex operators obey the operator algebra, similar to 
the operator algebra of the one-dimensional case (Sec. 3.4.1] ) . The difference is that the operator 
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algebra is localized in each point of the Fermi surface and the holomorphic coordinate is chosen to 
be relative to the Fermi momentum: 

V s (k f ,z)c a (z') ~ ( Z * f ~ Zkf r/ 2 c a (z')V s (k f ,z), (4.45) 

V c (k f ,z)c a (z') ~ { ^~ L Z ^ ) l / 2 c a {z')V c {k f ,z). (4.46) 

Below we also use a two-particle vertex operator [a(u, w) which creates a soliton and antisoliton of 
the spin displacement 

f x(u,w) = V s (u)V s - 1 (w). (4.47) 



{v! -uYl 2 {w' -wfl 2 



The operator algebra gives 

n(u', w')c\{u)c\{w) = — — — 

1 [u — W) 

x c\(u)c\(w)fi(u,w') (4.48) 

and if the points u, v! and w, w' coincide 

c-\(u)ci(w)fi(u,w) ~ _ n(u, w)c T (tt)c;(tD). (4.49) 

Two composite objects V s (z)c^{z) and V s (z)c[(j&) are singlets but carry electric charge. 

4.4 Electronic spectral function 

The bosonization and operator algebra provides a recipe for computing the most interesting Green 
functions. The calculations are similar to the one dimensional case (Sec. 3.3.2). 



We start from the spectral function. According to ( 4.36 ) the electronic Green function at equal 
time is 

n(r) = (4(r) C(J (0)} ~ J e ikr (k|y c t(k, r)K(k, 0)|k)(k|^(k, r)V a (k, 0)|k)dk, (4.50) 

where the state |k) = \l/(k)|0) is a zero mode state with momentum k. The first factor describes 
the amplitude of soft modes of density modulation in the state |k) whereas the second n-o(k) is 
the number of gapped excitations with spin 1/2. The spin excitations are gapped and their matrix 
elements are regular. At large r 3> vj//\q one can approximate the spin factor in fl4,50| ) by the 
occupation number of massive particles 

„ o(k )4(i- 7 =a£=M= ) . (4.5D 



so that the tunneling occupation number becomes 

n(r) = (4 (r) c CT (0)> ~ j e ikr (k|V;t(k, r)V c {k, 0)|k)n (A;), (4.52) 
The operator algebra gives the conformal block for the matrix element of the soft modes: 

1 e *6(k,r) 

Ac(r) = (k|K f (k,r)y c (k,0)|k) ~ — - = — — , (4.53) 
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where 9(k, r) is the angle between vectors k and r. As the result we have a "tomographic" repre- 
sentation for the occupation number 



n(r) ~ 

r 



In momentum representation the tunneling occupation number is 

n{k) ~ y D k (q)n (£ - 9 )dq, (4.55) 

where the propagator 

D k (q) = — = I e *(Argk-Argq) ,4 gg) 

q • k + ic[ x k q 
is a holomorphic function relative to k. 

As is in to the one dimensional case, the Green function may be obtained by replacing no(k) 
and -Dk(q) i n the eq. ( [4.55| ) by 

G (w,k) = a" 7^* (4.57) 

£>k(w,q) = 2 1 2 e <(Argk-A,gq) (4.58) 
w z — 

so that 

G(fi,k) ~ J G (n-OJ,k-Cl)D k (u,fl)dq,du. (4.59) 

and the spectral function 

A(k,E)~I ^(Argk-Argq) ^(k-q) 

•/ Uv}{k f + k - q)2 + A 2 - Ef - vjq 2 



where the integral goes over the domain E > vjq + yv 2 (k^ + k — q) 2 + A 2 , and the vector kj of 
the length kt is directed along k. 

The behaviour of the tunneling occupation number and the density of states close to the Fermi 
surface is linear ( modulo logarithmic corrections) 

k — kt 

n(k) -n(kf) ~ -const \k/k f - 1| sgn(/c - k f )0(\n— — £); at Vf(k - kf) » A . (4.61) 

k f 

As in the one dimensional case ( but contrary to the free massive particle case), the tunneling 
occupation number does not vanish at the Fermi surface. 

The tunneling density of states has a similar behaviour. It shows an asymmetric broad (~ e/) 
peak and decays from the peak toward the threshold E = Ao as 

dN(E) f _dft E E 



dE 



/JL Tp Tp 

A(k,E)^— ~ | — |C(ln — ), at e f » E - A » A (4.62) 
2ir €f €f 



In contrast to free massive particles but similar to the one dimensional case, the tunneling 
density of states has lost its singularity (E 2 — A 2 ,) 1 / 2 . It remains smooth at the threshold and 
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approaches it linearly (~ E — Ao). As in one dimension the spectral function is determined by 
two scales, Ao and e/, rather than just Ao- The second scale is the signature of the orthogonality 
catastrophe and of the composite nature of the electron. 

These formulas are to be compared with the one dimensional spectral function of the Sec, 3.3.2[ 



They are similar except that in one dimension the integral over the Fermi surface must be replaced 
by the sum over two Fermi points. While integrating over the Fermi surface special attention 
should be paid to the phase factor. This factor is not important in one dimensional case because 
the Fermi surface is not connected. In higher dimensions this factor establishes phase coherence 
between different points of the Fermi surface. 

4.5 Matrix elements 

The technique to compute the matrix elements also has been illustrated on the one dimensional 
case example (Sec. p.5.2 ). We adopt similar arguments for the two dimensional case to compute the 



two particle wave function: 

A(x-y)=e aa ,(N\c a (x)c CT ,(y)\N + 2). (4.63) 

4.5.1 Pair wave function 

First we must understand the nature of the ground state with two extra particles \M + 2). This 
state is created by a composite operator which creates a flux and places particles into zero modes 
created by the flux. Creating a flux also generates spin waves. They are gapped. Let us first 
assume that the gap Ao is very large and so that we can neglect spin excitations. Then the vertex 
operator V c of the charge sector creates a flux and a zero mode which gets occupied by a particle. 

Let us start from the two-particle state with zero momentum: 

|AA + 2)~ J (irdr / dke lk ( r - r 'V- 1 (k,r)F- 1 (-k,r / )^ t (k)^t(_k)|AA). (4.64) 

A physical meaning of the composite operator which creates the ground state with an extra 
particle is more transparent in terms of gauge invariant electronic operator ( 4,36j ). Using the 



tomographic representation we construct the two-particle state ( |4,64[ ) which is composed of the two 
electrons and the operator ( [4.47| ), creating a vortex and an antivortex of spin density at the points 
of electron insertions: 



\N + 2) ~ <w / drdr'dk (4.65) 

x V s (k,r)V7 1 (-k,r')4(k,r)( V (-k,r)|JV). 

Electrons with opposite spins see each other with vortices of the opposite angular momenta I = ±1 
(as it has been first noticed in Ref. p5| ). This gives an angular momentum / = 2 for the pair. 

The singlet two-particle matrix element has the form 

A(r) ~ J dke- 2iArg We ikr 
x(N + 2\V s (r)V s - 1 (0)\N + 2). 
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Thus the two-particle matrix element is given by the correlation function of the vertex operators 
of the spin channel computed in the ground state with two extra particles. The operator algebra 
V s (z)V~ l (z') ~ ^ : V s {z)V- l {z') : implies (J\f + 2| V a (k, r)^" -1 (k, 0)Af + 2) ~ ^e ie ^ r \ where 
9 (is., r) is the angle between vectors k and r. As a result 

A(r) ~ J_ e -2iArg(r) f g-i^k.r^ikr^ u . 66 ) 

kfr J 

The matrix element ( |4.66 ) is obtained in the limit of a very large gap Ao — ► oo. This approxi- 
mation is sufficient in order to analyze the transformation properties and the angular dependence of 
the tunneling amplitude. If the gap is not very large, an embedding of two extra particles creates a 
gapped spinon-antispinon excitations. These excitations do not interact with the zero mode of the 
charged sector and their wave functions (at kr 3> 1) are just the wave functions of an unperturbed 
theory flOg ) 

d>+(k, r) = u fc _ fc/ e* kr + v k . kf e^~ 2 ^\ (4.67) 
d> (k,r) = v k . kf e^ r - u k . kf e^- 2 ^>. (4.68) 

Here the first (second) function corresponds to a positive (negative) energy. Notice, that they are 
the same as in one dimensional case 

The wave function of a spinon - antispinon pair with opposite momenta is 

A BC s(k, n - r 2 ) = sign(A: - k f ) $+(k, ri)$„ *(k, r 2 ) (4.69) 



where the only translational invariant term retained. 

Let us notice that in contrast to Cooper's pairing mechanism, the gap in the topological mech- 
anism of superconductivity is generated via backward scattering and is of insulator nature. Never- 
theless, the wave function of a singlet spin excitation with zero relative momentum is the same as 
the BCS wave function, and is mainly determined by the Lorentz invariance. 

To obtain matrix elements for a finite gap, we have to replace the factor e ikr in eqs.(fh66) with 
the spinon-antispinon wave function (|4.69|) 



A(r) ~ e ~2iArg(r) 1 f e -ifl(k,r) A BCS (k, r)dk. (4.70) 
kfr J 

It is rather straightforward to take into account the frequency dependence of the two-particle 
matrix element. In order to do it one must replace -Dk(q) in the eq. ( 4.74j) by ( 4.5$ ) and Abcs by 

ABCs( "' k) = ^ _ (e(k) A _I! ^ _ ^ (4-71) 

so that 

A(fi,k) ~ e ~ 2iArg(k) J dqduj 

x A B cs(n-o;,k-q)D k (a;,q). (4.72) 
These results lead to a number of important consequences. 
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4.5.2 Orthogonality catastrophe and angle dependence of the tunneling amplitude. 
Complex d-wave state 



(i) The angular dependence of the pair wave function and of the tunneling amplitude. The pair 
wave function consists of two vortices of the same charge — one comes from the wave function of the 
zero mode of the flux created for a spin up particle by a spin down particle (the factor e _jArg ( r )^ 
in ( 4.70D ). Another vortex (e~* Arg ( k )) is located at the center of the Fermi surface. The pair 



wave function forms the d-wave (I = 2) irreducible complex representation of the rotational group. 
Similarly the one particle matrix element carries / = 1 orbital moment. The angular dependence 
of tunneling in anyon superconductors has been also studied in Refs.||]. 

(ii) Orthogonality catastrophe. The overlap between two ground state wave functions with M 
and J\f + 2 electrons does not vanish in a large system. 

This is not the case for a one- (or any odd number) particle matrix element. An attempt to 
embed a single electron in the system leads to a half-occupied zero mode state. As a result this state 
turns out to be almost orthogonal to all other states with the same spin and number of particles. 
The matrix element (M\c a (k + q)\M + 1, k, q) vanishes at q — » 0. 

(Hi) Tomographic representation. Eq. ( [4.70 ) consists of the integral over the Fermi surface and 



may be viewed as a "tomographic" representation of the matrix element 

A(r) ~ e ~ i2A ^ J dkD k (r)A BC s(k, r) (4.73) 

where the propagator -Dk(r) is given by ( [4.53| ). The tomographic representation in electronic liquids 
has been anticipated in p6| , p7| ]. The new feature is that the relative phase of electron pairs at 
different points of the Fermi surface are correlated by the factor e lArgk in the propagator of soft 
modes (|4.53|) . 



(iv)Bremsstrahlung. It is instructive to rewrite eq. ( 4.73 ) in the momentum representation 



A(k) ~ e - 2A r g ( k ) J A BC s(k - q)Ac(q)dq (4-74) 



where the propagator £*k(q) is given by (|4.56|) and 



A BC s(k) = . A ° (4.75) 
/(e(k) - /i) 2 + A 2 



Similar formula stays for the tunneling occupation number (4.55). 

These representations clarify the physics of the topological superconductor. An insertion of 
two particles in the spin singlet state with relative momentum k close to k^ emits soft modes 
of transversal spin current with the propagator -Dk-kj-(q)- As a result of this: (i) ground states 
differing by an odd number of particles are orthogonal (for a related arguments for a one dimensional 
spin chain, see p3j); (ii) the BCS wave function is dressed by soft transversal modes. This is in 
line with ID physics and the bremsstrahlung effect of QED. The new feature is the phase of the 



matrix element of the soft mode (4.56) which is the angle relative to the Fermi momentum of the 
electron. By contrast, in the BCS the density modulations are not individual excitations but are 
composed from Cooper pairs. The interaction between density modulations and pairs — the major 
effect of the topological mechanism — vanishes in the BCS. 
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(iv) Momentum dependence of the two-particle wave function and tunneling amplitude. The 
momentum dependence of the amplitude of the pair wave function |A(k)| is drastically different 



from the BCS (4.69). In the vicinity of Fermi surface \k — kf \ <C kf the integrals (4.74) can be 
computed: 

\A(k)\ -\A(kA\ ~ i v f( k ~ k f)/ A o v f \k-kf\ < A 

1 1 jl 1 { f)l ~ \ sgn (jfe - k f ) log(v f \k - fe/|/Ao) v f \k - k f \ > A 

The result gives the universal dependence of the pair wave function on Vf(k — kf) /Aq. The constant 
|A(fcf)| is not universal and depends on states far away from the Fermi surface. 

In contrast to the BCS gap function, the pair wave function is asymmetric around the Fermi 
surface. It is peaked at scale 

Vf(k- kf) « €f ( lQ g^) » 

which is much greater than Ao- Also, the pair wave function has a logarithmic branch cut at 
Vf(k — kf) = HAq in contrast to square root singularity of the BCS function. This indicates that 
transversal spin current soft modes are emitted in the process of tunneling. 

4.6 Pair Correlation Function and Long Range Order 

A non vanishing two particle wave function does not mean that the pair correlation function 

K(n,r2,r 3 ,r 4 ) = (cj(ri)cj(r 2 )c t (r 3 )c|(r4)) (4.76) 
shows long range order. We already observed this result in one dimensional case (Sec. (13X11) ). 



Indeed, if we assume that size of a pair is less than the distance between pairs |rx — r 2 | ~ 
|r 3 — <C |R| = |ri + r 2 — r 3 — and consider the dependence of the correlation function on R, 
we have 

K = /(c{(k,ri)c{(-k,r 2 )c T (p,r 3 )c i ((-p,r4)))dkdp ~ 

J dke ik ( ri ~ r2 )e~ ik ( r3-r4 ) ( e -*c(k,r 1 )-* c (-k,r 2 )+*c(k,r 3 )+* c (-k,r4)^ x 

/ e -$ s (k,r 1 )+# s (-k,r 2 )+# s (k,r 3 )-0 s (-k,r4)\ (4.77) 

As in the one dimensional case, the spin sector factor approaches a constant at R — > oo, whereas 
the charge sector factor decays. As a result the pair correlation function decays as well: 

K(r 1 ,r 2 ,r 3 ,r 4 ) ~ A / e ik ( ri ~ r2 )e^ ik ( r3 - r4 )e 2i ^ k ' R )(ik (4.78) 
R 2 J 

In contrast to the BCS, the decay of the pair correlation function does not mean that a supercon- 



ducting state is destroyed by fluctuations. Nor does it mean that the tunneling amplitude (4.63) 



vanishes (see the comment in the end of Sec. 3.3.4 ). The decay rather indicates a strong interference 



effect between pairs. Similar to the Berezinsky-Kosterlitz-Thouless transition, the rate of the decay 
reflects the value of the superfluid density. 

It is easy to construct an object which carries the electric charge 2e and shows a long range 
order. As in one dimensional case (see ( |3,74| )) it can be written as ^(r\,r2) = 
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/ c{(k, ri )V c (k, ri )V c (k, r 2 )c{(k, r 2 )dk 

K(R) = (*t( ri ,r 2 )*(r 3) r4)) ~ 

J e ik(ri-r2) e -ik(r3-r4)^ e -$ s (k,r 1 )+* s (-k,r 2 )+*c(k,r 3 )-$ s (-k,r 4 )^ ( ^ c ^ i &t R 



OO. 



(4.79) 



This long range order indicates a gap opening at 2k j rather than a superfluid order. 

A striking difference to the BCS is a relation between the phase of the tunneling amplitude 
and the phase which parameterizes the current ( |2.19| ). In the BCS these phases are the same. In 
topological superconductor there is no direct relation between the phase of two particle matrix 
element ( |4,63 ), the pair correlation functions ( 4.76| ) and ( [4. 79 ). 



5 Tunneling and orthogonality catastrophe 

Tunneling effects are most peculiar in the presence of an orthogonality catastrophe. Tunneling 
experiments seem to be the best instruments to search for a signature of the topological mechanism. 
Below we briefly discuss various tunneling effects. 

5.1 Josephson Tunneling 

The Josephson current through the junction (at zero bias voltage) between two superconductors, 
of which one is conventional, is given by the well known formula JT3| 

x , A? 

I = -Smi^TkpT-k-p— 

kp k 



00 du 
2tt 



F(p,u) F(p,-u) 



E k -ir] to- E k +irj 



(5.1) 



where T kp is a transmission amplitude of the junction, E k = J | Ai| 2 + is the pair spectrum of the 
conventional superconductor and F(p,uj) is the spectral function of the superconductor of interest 

F(p, w) = 27ri£ CTCT ' ^2(N\c a (p)\N + 1){/V + 1|<v(-p)PV + 2)<J(o; - e p ). (5.2) 

Here the sum goes over all quantum states with one extra particle (e p is the energy of an intermediate 
state). If the spectrum is symmetric with respect to adding or removing a particle, i.e., F(p,co) is 
an odd function of a;, we obtain 

I = -2<5mi^T kp T_ k _ p — / V- —■ (5.3) 

^ y ' E k Jo 2vr lo - E k 

In the following we assume for simplicity that the transmission amplitude T kp is strongly peaked 
at k, p close to the direction normal to the junction. This simplification should not change 
the phase dependence of the Josephson current although it can change the value of the critical 
current. Assuming that the gap in the superconductor of interest Ao is bigger than the one in the 
conventional superconductor Ao S> Ai we obtain: 

poo 

/~ llfz/olAxI V / — sin((fo-^(p n ,u;))|F(p n ,u;)|, (5-4) 
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where p n is the component of momentum normal to the surface of the junction and averaging over 
p n is determined by the actual form of the transition amplitude. In this formula 4>q and vq are 
the phase and density of states of the conventional superconductor and </>(p, u>) is the phase of the 
F- function (5.2) of the superconductor of interest. 



In the BCS theory the F- function has a peak at the gap u ~ Ao, such that the width of the 
peak is also of the order of the gap: 

dpF B cs(p,u) r A ° (5.5) 



The peak gives the major contribution to the integral ( |5.4| ). It selects a characteristic energy of 
the intermediate state e q ~ Ao and gives rise to the traditional BCS picture of tunneling: a pair 
decays into two electrons while tunneling, so electrons tunnel independently. Short time processes 
(at co ~ ej) do not contribute to the integral 



The situation is drastically different in the orthogonality catastrophe environment [11|. In a 
topological superconductor an individual matrix element (j\f\c a (p)\j\f + 1) acquires an additional 
factor 1/Lkf, where L is the size of the system and therefore vanishes in a macroscopical sample. 
In other words the ground states with j\f and j\f + 1 particles are almost orthogonal. Nevertheless, 
the tunneling, i.e., a matrix element between states with j\f and N + 2 particles is nonzero due to 
a large number of low-energy intermediate states contributing to the sum fl5.2j). A result of this 
is that the pair spectral function (|5.2| ) acquires an additional factor uj/tf- In contrast to BCS at 
uj S> Ao we have 

/ dpFM ~ ^ (-) O(ln^) ~ ^O(ln^). (5.6) 
J uj \ej J A e/ A 

Therefore, the characteristic scale of the spectral function is shifted to the ultraviolet and becomes 
of the order of the Fermi energy — much larger than the scale of the gap : the integral ([5^4]) is 
saturated by u ~ ef ^> Aq. This means that a pair remains intact during the tunneling and the 
tunneling amplitude is determined by the equal time value of F(p, t = 0), i.e., by the matrix element 
A(r) of an instantaneous creation of a pair ( |2.26| )P| Let us notice that the correction to the spectral 
function in the eq. ( |5.6[ ) — {oj/ej) a at a = 1 is just marginal. Were a is less than 1 the time of the 
tunneling would be of the order of Aq 1 . 

In the coordinate representation the Josephson current is 

/ ~ 9m e"*A(fe/n) ~ sin(0 o - <j>(k f n, u = 0)). (5.7) 

In the corner-SQUID-junction geometry |9j one can directly measure the difference of phases 
of A between two faces of the superconducting crystal. If the faces of two contacts have a relative 



angle 9, then, according to ( 4.73 , 4.74 ) the phase shift in tunneling amplitude will be twice the angle 
0(fc/ni) - <t>{k f n 2 ) = 26: 

A ( fc / n i) _ p 2ie ( , R) 

ATv^y • ( } 

The instantaneous character of tunneling and the power laws ( |5.6| ) are known from one-dimensional 
electronic systems where the orthogonality catastrophe comes to its own. In addition, in 2D it also 
leads to the angular dependence of the two-particle amplitude. 

8 A similar phenomenon has been discussed by Chakravarty and Anderson in the context of interlayer tunneling 
in cuprate superconductors 
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5.2 One particle tunneling 



The one particle tunneling current is also strongly affected by the orthogonality catastrophe. One 
particle tunneling current from a metal or a superconductor with a small gap is determined by the 



tunneling spectral function (4.62) 



, ,o f eV dEdN(E) , . 

/ "' = 2|T| "»A ^-hr* (5 - 9> 

where vq is the density of states of the normal metal (or low gap superconductor). In the BCS 
dN J^^ ~ ~UW A 2 anC ^ ^ e ^^ rec ^ curren t is given by a familiar formula I ~ \J (eV) 2 — A 2 ,. 

In the topological superconductor the one-particle Green function acquires a branch cut, rather 
than a pole on the threshold and the tunneling density of state looses its singularity at the threshold 
( |3.62j ). This leads to a suppression of direct current close to the threshold bias voltage. With a 
logarithmic accuracy we will obtain almost quadratic I — V behavior well above the threshold 

~ V, at eV > A . (5.10) 

It peaks at Aolne//Ao- 
5.3 N-S Tunneling 

The conductance between normal metal and a superconductor does not vanish for the voltage 
below the gap due to the pair tunneling: an electron incident on the NS boundary captures another 
electron and generates a hole in the metal and pair in the superconductor (Andreev scattering). If 
one assumes that tunneling Hamiltonian for the pair tunneling is 

H NS = rC T (r 1 )q(r 2 )4(r 3 )c{(r 4 ), 

than the tunneling current of this process is given by 

I = e\T\ 2 v$ J V duj{eV - u>) J QmK(uj, q)dq (5.11) 

where 

K(u,q) = J e i ^ +iqR (c{(ri,t)c{(r + ri,t)c T (R + r 1 ,0)q(R + r / + ri,0))(iR^ (5.12) 

is the pair correlation function. Its dependence on positions of electrons r, r' within pairs is taken 
care by the tunneling matrix element. At small voltage only small pairs contribute r,r' <C R to the 
tunneling current. 

In the BCS, the pair correlation exhibits a long range order, i.e., QmK(ui,q) ~ p s 5(q}5(uj), 
the integral in ( |5.11 ) is saturated at q, to = and gives a nonzero conductance I ~ &qV, where 



a ~ (|r| 2 e 2 z/ 2 p s ). 

In contrast, in the topological superconductor, the pair propagator is determined by soft col- 
lective density fluctuations. Extending the result of Sec. 4.6, we get K(u, q) ~ (uj 2 — v 2 q 2 ) 1 ^ 2 . This 
gives a nonlinear I — V characteristic 

eV 

I ~ a V( — ), at eV <C A . 

e f 
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This result implies zero NS conductance and a suppression of the Andreev scattering by the or- 
thogonality catastrophe 



6 Concluding Remarks 

6.1 Physics of one dimension beyond the one dimension. 

The goal of this paper is to show that one dimensional electronic physics is not restricted to one 
spatial dimension. The origin of almost all phenomena we have observed in one dimension is the 
spectral flow, rather than a restricted geometry. A very similar physics takes place in topological 
electronic liquids, i.e. in compressible liquids with a charge spectral flow. 

One of the features of one dimensional physics - the Sugawara form of the stress energy tensor 
and the associated current algebra, being lifted to a higher spatial dimension, translates into su- 
perconductivity. With a great deal of generality we argued that a spectral flow in a compressible 
liquid with a fixed chemical potential inevitably leads to the hydrodynamics of an ideal liquid ( 
a generalization of the Sugawara form of one dimension) and to a superconducting ground state 
(topological mechanism of superconductivity) . 

Another lesson may be learned from one dimension - the orthogonality catastrophe is also a 
consequence of the spectral flow. It affects drastically all processes of the topological superconductor 
involving a change of number of particles. Among them one particle spectral properties (e.g. 
photoemission) and all sorts of tunneling (Josephson tunneling, NS two particle tunneling and 
Andreev scattering, one particle tunneling). 

Neither electrons nor Cooper pairs are elementary excitations in topological liquids. The in- 
sertion of one electron (any non-singlet state) drastically changes the ground state of the system, 
so that the matrix element between two ground states with N and N + 1 electrons vanishes in 
a macroscopic system. The same is true for any ground states which differ by an odd number of 
particles. On the contrary, the matrix elements of two particles in a singlet state do not vanish but 
are significantly modified by the interaction with soft modes of density modulations. In particular, 
the poles of Green functions on the mass shell are replaced by branch cuts. 

The one-dimensional physics provides tools to compute matrix elements and Green functions 
in topological fluids in dimension greater than one. This is bosonization. Actual particles, having 
asymptotic states, are solitons of nonlinear charge density and spin density modulations. They 
can be represented as coherent states of the charge and spin densities. Their matrix elements are 
composition of holomorphic functions, no matter what the spatial dimension is, and can be treated 
in a similar way to what used in one dimension. One is the method is bosonization, presented in 



Sec. 4.3. Another (more general) method is to consider the matrix elements as conformal blocks 
of a certain conformal field theory and a current algebra. Let us notice that electronic liquids 
without a spectral flow, say the Fermi liquid, can not be bosonized, i.e. their bosonized form is not 
adequate to physics of Fermi liquid ground state. The reason for this is a considerable phase space 
of scattering at the Fermi surface (leading to a dissipation). Contrary in topological liquids the 
scattering and emission of a soft mode are asymptotically forward. This reduces the phase space 



9 This result will not be affected by the Coulomb interaction, provided that plasmon frequency is smaller than 
the gap. Effects of interaction between Cooper pairs and soft modes of density modulation on various tunneling 
mechanisms in conventional BCS and layered superconductors has been studied in 1291] . 
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of scattering and eliminates dissipation. On a general basis one can say that only superconductors 
can be bosonized. 

6.2 Topological superconductor versus BCS superconductor 

The orthogonality catastrophe is a source of major differences between a topological and BCS 
superconductors: 

(i) The tunneling amplitude or the ground state pair wave function A(ri, Y2) = (AA+2|cj(ri)cJ(r2)|A/ r ) 
differs from the gap function. The former is complex and its phase 2Arg(ri — depends on the 
direction of the pair, whereas the gap may be isotropic. The amplitude of the pair wave function 
and the gap function are also different. The scale of the former is determined by the gap and also 

by the Fermi scale. 

(ii) Contrary to BCS the gradient of the phase of the pair wave function in a current state 
A( r i) r 2) = |A|e l<,p does not determine the current j 7^ — ^—^7ip. The relation between the phase 
and the current is nonlocal and retarded. The pair wave function does not play the role of an order 
parameter. 

(hi) The pair correlation function (cj(ri)cj(r + ri)c-f (R + ri)cj_(i? + r' + ri)) does not have a 
long range order but decays with the separation between pairs at R — > 00. Nevertheless, the 
superffuid density and the tunneling amplitude do not vanish. 

(iv) The pair spectral function A(k) has a broad structure of the amplitude of the pair wave 
function in momentum space around what one may call Fermi surface. The width of this structure 
is of the order of fc//log(e//Ao) which is much bigger than Aq/vj, the width of the peak of the 
BCS wave function. 

(v) The origin of the gap is the 2k f instability rather than a condensation of pairs as in BCS. 

(vi) Tunneling: The time of Josephson tunneling is determined not only by the gap (as in BCS) 
but also by the Fermi scale and is much shorter than Ag 1 . A pair remains intact while tunneling. 

One particle tunneling is suppressed above the threshold. The amplitude of Andreev scattering 
vanishes at the Fermi surface, so that a conductance of the NS tunneling. /- V characteristic of the 
NS tunneling becomes nonlinear. 

(vii) The noise spectrum (j(u;)j(— uj)) of a point Josephson junction is expected to have a power 
law peak, rather than a narrow peak of the BCS superconductors (we do not study the noise 
spectrum here). 

(viii) One particle spectral function has a branch cut rather than a pole i.e., Z-factor (a residue) 
as well as the one particle matrix element vanishes at the Fermi surface. The one particle spectral 
function (tunneling density of state) has a broad peak shifted inside the Fermi surface and vanishes 
with a (weak) singularity at the Fermi surface. 

Ironically, the quantities listed above did not enjoy a detailed measurements even in standard 
superconductors. 

6.3 Parity and time reversal symmetry breaking. 

There have been misconceptions in the literature regarding parity symmetry breaking of the ground 
state of two-dimensional topological liquid. The spatial parity and time reversal symmetry are 
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simultaneously broken in the ground state of the model considered. This reflects the chiral nature 
of zero modes. However the broken symmetry does not exhibit itself in every physical quantity. 
Also there is no easy ways to detect this symmetry breaking experimentally, if any at all. The 
reason is that, although the time reversal symmetry is broken, there are no spontaneous local 
electric currents either in the bulk or on the edge in any steady state. One can see this from the 
hydrodynamics of the charge sector ( 2.16| ), but this fact remains valid beyond the low frequency 
range. Even more general, all diagonal singlet matrix elements are parity and time reversal even. 
Therefore, one should not expect to observe the time reversal symmetry breaking by measuring 
Faraday rotation |^] and muon spin relaxation [31] — early experimental searches for a signature 
of parity symmetry breaking. 

Off-diagonal or non-singlet matrix elements are a different matter. The broken time reversal 
symmetry is explicit in the complex d-wave tunneling amplitude. Another manifestations of the 
broken symmetry can be seen in the spin sector. Among them are an expectation value of the spin 
chirality S • VS x VS and a novel feature — edge spin current. This is a two-dimensional version of 
the known phenomenon in ID. A spin chain with gapped bulk spin excitations develops gapless spin 
excitations at the edges. In a two-dimensional spin liquid edge excitations are chiral spin currents. 
Edge magnetic excitations have been observed in spin chains ]32[ |. One may expect to find these 
soft edge spin excitations in model systems with an enhanced boundary (say, for example, an array 
of superconducting islands). 

Parity symmetry breaking is not an inherent property of the topological fluids, but rather a 
property of the "irreducible" model. By doubling the number of fermionic components one can 
easily construct a model with no time reversal or parity symmetry breaking. The most natural 
way to do this is to include an even number of layers. The parity alternates between odd and even 



layers p3, 18 1. Genuine parity breaking takes place only in systems with odd number of layers and 



minimal number of fields. 



6.4 Topological mechanism and superconductivity of cuprates 

Being a theory of doped Mott insulators, topological superconductivity has a number of features 
observed in cuprate superconductors. A detailed discussion of the experimental data from the point 
of view of topological superconductor is far beyond the scope of this text. We only mention the 
corner-SQUID-junction experiment Q in which the relative phase of tunneling amplitudes (2.23) 
on different faces of a single crystal of YBCO has been measured. It is found to be tt in accord with 
d-wave superconductivity and apparently in agreement with the topological mechanism ( 4.1 9| ) . 

A critically inclined reader faces difficulties to determine weather known London type super- 
conductors are BCS superconductors or topological liquids. 
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7 Appendix. Fermionic number and Solitons 



For completeness, let us review a well known calculation of the fermionic number (3.4, |3.16| ). It is 



(also known as Atiah-Patodi-Singer invariant 77 ) is a number of unoccupied states appeared as the 
result of adiabatically created topological configuration of the potential. 

Since without a potential the spectrum to the Dirac operator is symmetric, this number also 
measures the difference between the number of states with positive and negative energy ij = A/+ — 
AA_ (spectral asymmetry). If AA = AA + + AA_ is the total number of levels, then the number of 
negative levels would be AC = (AA — tj)/2. The "regularized" number of negative levels, i.e., the 

2 2' 



change of the negative levels due to the potential, AAA = AA_ — ^r = — 1?7, is the fermionic number: 



AAA = -^Trfsign H] . (7.1) 

This divergent quantity has to be regularized. One of the standard ways of regularization is to 
replace it by 

The integrand now converges. 
One Dimension. 

We compute the fermionic number for a more general Peierls model 

H = a x id x + 0m + ia x /3TT 2 (7.3) 

where the Dirac matrices a x , (5 may chosen as the Pauli matrices: a x = 03, = a\. We also assume 
that the modulus of the vector (tti,^) takes a fixed value at infinity. 

The soliton here is a field n(x) which forms a homotopy class 7Ti(S 1 ). Its topological charge is 

(7.4) 

With a given value of the |7r| at infinity Q is an integer number. 
The square of the Hamiltonian is 

H 2 = (id,) 2 + vr 2 + e ab a a d x 7r b , a, b = 1, 2 . (7.5) 



Let us now expand the integrand of ( |7.2| ) in dir. In the first non- vanishing order of |7?(oo)| — > 00 
we have 

1 f +co _ ( a 3 id + aTT . 3 1 



AAA = — — / Tr — — ^ia-'adn-—, ? (7.6) 

2^7-00 \(id) 2 + IT 2 + z 2 (id) 2 +vr 2 + z 2 ) v ; 

This is the only order which contributes to the fermion number. The trace in spinor space gives 

1 r + °° , nn ( ^dir 2 - ^dir 1 \ 



The trace in momentum space is conveniently calculated in the plane wave basis, 

I dx I d i e ^WT^TW ' (7 ' 8> 
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Finally the integrals over p and z give the topological charge ( |7.4|) of kinks 

AM = -—Jdx e^—^- = Q , (7.9) 

the topological charge ( |7.4[) of kinks. The number of extra states induced by a soliton therefore 
equals the topological charge of the soliton. If tt\ + iiri = Aoe l</3 , we obtain a global version of the 
Frohlich relation ( |3.12[) . The model with 7Ti = deserves special interest. This is commensurate 
Peierls model (iT2 = A). It is defined by 

H = a x id x + i l5 A . (7.10) 

The spectrum of the Hamiltonian is symmetric. This means that if there is a state ipE with an 
energy E then there is always a state (3ip with the energy —E, except for E ^ 0. It follows from 
the fact that the Hamiltonian anticommute with the matrix (5. 

Now the soliton is a kink with A(— oo) = — A(+oo) = Ao- Setting iti = in Q7.9| ), we obtain 

Q = -^J 5(A)dAdx . (7.11) 

In this case an extra state appears in the middle of the gap (zero mode). The kink does not respect 
periodic boundary conditions. As a result of this a zero mode state may accommodate only 1/2 
of the particle. In a system with a periodic boundary conditions kinks may appear only together 
with antikinks, so the total number of states remains an integer. 



Two Dimensions. 

Calculations are similar in two dimensional case 

H= 4[a(iV + A)+/3n# ff . (7.12) 

(7=1,2 

The square of its Hamiltonian is 

H 2 = (iV + A) 2 + m 2 - /3(V x A). (7.13) 



Let us expand the integrand in ( |7.2j ) in terms of /3(V x A). The first and only non- vanishing order 
gives: 

a(iV + A) + f3m 

.(iV + A) 2 + m 2 + z 1 

1 



AiV = -±/ + ^Tr J (VxA)x 



(iV + A) 2 + m 2 + z< 

The trace leaves the integral over p and z: 



(7.14) 



2 f + , f ,2 f d P m(V x A) _ . 

AN = / dz / d 2 X / ^-Wtt 7T oTn ■ (7.15) 

ttJ-oo J J (2vr) 2 [(p + A) 2 + m 2 + z 2 ] 2 V ; 

As a result we obtain the relation between the fermion number and the topological charge of the 
vortices (fl2|, |4~3|) 

AiV = sign m / d 2 x (V x A) . (7.16) 

2-7T J 
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